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Abstract 

We show that the apphcation of the pinch technique to the conventional Schwinger-Dyson equa- 
tions for the gluon propagator, gluon-quark vertex, and three-gluon vertex, gives rise to new 
equations endowed with special properties. The new series coincides with the one obtained in the 
Feynman gauge of the background field method, thus capturing the extensive gauge cancellations 
implemented by the pinch technique at the level of individual Green's functions. Its building blocks 
are the fully dressed pinch technique Green's functions obeying Abelian all-order Ward identities 
instead of the Slavnov- Taylor identites satisfied by their conventional counterparts. As a result, 
and contrary to the standard case, the new equation for the gluon self-energy can be truncated 
gauge invariantly at any order in the dressed loop expansion. The construction is streamlined 
by resorting to the Batalin-Vilkovisky formalism which allows for a concise treatment of all the 
quantities appearing in the intermediate steps. The theoretical and phenomenological implications 
of this novel non-perturbative framework are discussed in detail. 
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I. INTRODUCTION 

The Schwinger-Dyson equations (SDE) provide a formal framework for tackling 

physics problems requiring a non-perturbative treatment. The SDE constitute an infinite 
system of coupled non-linear integral equations for all Green's functions of the theory, and 
can be used, at least in principle, to address questions related to chiral symmetry breaking, 
dynamical mass generation, formation of bound states, and other non-perturbative effects 
3|, |j| . In practice, their usefulness hinges crucially on one's ability to devise a self-consistent 
truncation scheme that would select a tractable subset of these equations, without compro- 
mising the physics one hopes to describe. Inventing such a scheme for the SDE of gauge 
theories is a highly non-trivial proposition. The problem originates from the fact that the 
SDEs are built out of unphysical Green's functions; thus, the extraction of reliable physical 
information depends critically on delicate all-order cancellations, which may be inadvertently 
distorted in the process of the truncation. Several of the issues related to the truncation of 
the SDE= of QED W been add.e.ed . a ot a.icle. flBB,! 

The situation becomes even more complicated for strongly coupled non-Abelian gauge 
theories, such as QCD 9|, mainly for the following reasons. 

i. The complications caused by the dependence of the Green's functions on the gauge- 
fixing parameter are more acute in non-Abelian gauge-theories, as can be seen already 
at the level of the most fundamental Green's function, namely the two-point func- 
tion (propagator) of the corresponding gauge bosons. In QED the photon self-energy 
(vacuum polarization) is independent of the gauge-fixing parameter, both perturba- 
tively (to all orders) and non-perturbatively; when multiplied by it forms a phys- 
ical observable, the QED effective charge. In contradistinction, the gluon self-energy 
is gauge-dependent already at one loop; depending on the gauge-fixing scheme em- 
ployed, this dependence may be more or less virulent. This difference is of little 
practical importance when computing ^-matrix elements at a fixed order in pertur- 
bation theory, given that the gauge- dependence of the gluon self-energy is guaranteed 
to cancel against similar contributions from other graphs, but has far-reaching conse- 
quences when attempting to truncate the corresponding SDEs, written in some gauge. 
Contrary to what happens in the perturbative calculation, even if one were to put 
together the non-perturbative expressions from these truncated SDEs to form a phys- 



ical observable, the gauge-cancellations may not go through completely, because the 
process of the truncation might have distorted them. Thus, there is a high proba- 
bility of ending up with a residual gauge-dependence infesting one's non-perturbative 
prediction for a physical observable. 

ii. In Abelian gauge theories the Green's functions satisfy naive Ward Identities (Wis): 
given the tree-level WI, the all-order generalization is obtained by simply replacing 
the tree-level expressions by the all-order ones. In general, this is not true for the 
Green's functions of non-Abelian gauge theories, where the Wis are modified beyond 
tree-level, and are replaced by more complicated expressions known as Slavnov- Taylor 



identities (STIs) [10|, lUl]; in addition to the basic Green's functions of the theory they 
involve various composite ghost operators. In order to appreciate how the fact that 
the Green's functions satisfy STIs may complicate the truncation procedure of the 
SDEs, let us consider the simplest STI (and WI in this case) satisfied by the photon 
and gluon self-energies alike, namely 

g"n,^(g) = 0. (1.1) 

Eq. (11.11) is without a doubt the most fundamental statement at the level of Green's 
functions that one can obtain from the BRST symmetry |l^; it affirms the transver- 
sality of the gauge-boson self-energy, be it a photon or a gluon, and is valid both 
perturbatively to all orders as well as non-perturbatively. The problem stems from 
the fact that in the SDE of Uapiq) enter higher order Green's functions, namely the 
fully-dressed fundamental vertices of the theory. It is these latter Green's functions 
that in the Abelian context satisfy Wis whereas in the non-Abelian context satisfy 
STIs. Thus, whereas in QED the validity of Eq. (11. ip can be easily seen at the level 
of the SDE, simply because gT^(p, p + q) = e [S^^{p + q) — S^^{p)], in QCD proving 
Eq. (II. ip is very difficult, and requires the conspiracy of all full vertices appearing in 
the SDE. Truncating the SDE naively usually amounts to leaving out some of these 
vertices, and, as a result, Eq. ( II. ip is compromised. 

The complications stemming from the two points points mentioned above are often com- 
pounded by additional problems related to the loss of multiplicative renormalizabihty and 
the inability to form renormalization-group invariant quantities. 
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Recent 
proposed 



V, a truncation scheme for the SDEs of non-Abehan gauge theories has been 



13] that is based on the pinch technique (PT) Q, and its connection with the 



background field method (BFM) |l6l.ll7l] (see below). The way the PT resolves the difficulties 
related with points (z) and (ii) mentioned above is by imposing a drastic modification already 
at the level of the building blocks of the SD series, namely the off-shell Green's functions 
themselves. The PT is a well-defined algorithm that exploits systematically the symmetries 
built into physical observables, such as S-matrix elements or Wilson loops, in order to 
construct new, effective Green's functions, endowed with very special properties. The basic 
observation, which essentially defines the PT, is that there exists a fundamental cancellation 
between sets of diagrams with different kinematic properties, such as self-energies, vertices, 
and boxes. This cancellation is driven by the underlying BRST symmetry [12], and is 
triggered when a very particular subset of the longitudinal momenta circulating inside vertex 
and box diagrams generate out of them (by "pinching" internal lines) propagator-like terms. 
The latter are reassigned to conventional self-energy graphs, in order to give rise to the 
aforementioned effective Green's functions. These new Green's functions are independent 
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of the gauge-fixing parameter 



of the complicated STIs 
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18|,ll9|,l20(, satisfy ghost-free, QED-li 



181 ] ■ display only physical thresho 
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i.e Wis instead 
, have correct 



2M- 



analyticity properties [23[, and are well-behaved at high energies 

Of central importance for the ensuing analysis is the connection between the PT and 
the BFM. The latter is a special gauge-fixing procedure that preserves the symmetry of 
the action under ordinary gauge transformations with respect to the background (classical) 
gauge field A^, while the quantum gauge fields, A^, appearing inside the loops, transform 
homogeneously under the gauge group [25|. As a result, the background n-point functions 
{i.e., those involving fields) satisfy QED-like Wis to all orders. The BFM gives rise to 
special Feynman rules (see Appendix [F|) ; most notably (a) the tree- level vertices involving 
A^j^ fields depend in general on the quantum gauge- fixing parameter ^q, used to gauge-fix the 
quantum gauge fields, and (6) the ghost sector is modified, containing symmetric gluon-ghost 
vertices, as well as two-gluon-two-ghost vertices. Notice an important point: the background 
ra-point functions are gauge-invariant, in the sense that they satisfy (by construction) QED- 
like Wis, but are not gauge-independent , i.e., they depend explicitly on Jn. The connection 



between PT and BFM 



26 



27|, demonstrated to be valid to all orders 



affirms that the 



(gauge-independent) PT n-point functions coincide with the BFM n-point functions when 
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the latter are computed at the special value = 1, also known as the background Feynman 
gauge (BFG). 

Let us now return to points (i) and (ii) and analyze them from the perspective of the PT. 
In a nutshell, the way point {i) is resolved, for the prototype case of the gluon self-energy, 
is the following. The BFG is a privileged gauge, in the sense that it is selected dynamically 
when the gluon self-energy is embedded into a physical observable (such as an on-shell test- 
amplitude). Specifically, the BFG captures the net propagator-like subamplitude emerging 
after QED-like conditions have been replicated inside the test-amplitude, by means of the 
PT procedure. Thus, once the PT rearrangements have taken place, the propagator is 
removed from the amplitude and is studied in isolation: one considers the SDK for the 
background gluon self-energy, Uapiq), at = 1. Solving the SDK in the BFG eliminates 
any gauge- related exchanges between the solutions obtained for Ilaf3{q) and other Green's 
functions, when put together to form observables; thus, the solutions are free of gauge 
artifacts. Regarding point (ii), the key ingredient is that now all full vertices appearing in 
the new SDK satisfy Abelian Wis; as a result, gluonic and ghost contributions are separately 
transverse, within each order in the "dressed-loop" expansion. Thus, it is much easier to 
devise truncation schemes that manifestly preserve the validity of Eq. (11.11) . 

The main results presented in this article are the following. We provide a detailed and 
complete demonstration of how the application of the PT algorithm at the level of the con- 
ventional SD series leads to a new, modified SD series, with the special properties mentioned 
above. A preliminary discussion of this issue has already appeared in brief communication, 
dedicated to the SDK of the gluon self-energy ISj. From the technical point of view, here 
we present a significantly more concise and direct proof, by virtue of a crucial STI, that 
we employ for the first time. In addition, we extend the analysis to include the SDK for 
the quark-gluon and three-gluon vertices, which are important ingredients for obtaining a 
self-contained picture. We emphasize that the three-gluon vertex relevant in this analysis is 
the one that would correspond, in the BFM language, to F^^^, i.e. one background gluon 
and two quantum ones merging, and not the fully Bose-symmetric vertex F^^^ considered 



m 
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29| . The reason is that it is the former vertex that appears in the SDK o: 



the gluon 



self-energy, in complete accordance with both the PT unitarity arguments 30| (see also 
Appendix and the independent diagrammatic rules of the BFM 17 1. 

Furthermore, we address an important conceptual and practical issue, related to the fact 
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that, qualitatively speaking, the new SD series expresses the BFG Green's functions in terms 
of integrals involving the conventional ones. This is already evident at the two-loop level: the 
two-loop BFM gluon self-energy is written in terms of integrals involving the conventional 
one- loop gluon self-energy. This example might suggest, at first sight, that one cannot arrive 
at a genuine SDE involving the same unknown quantity on both sides, but there is a way 
around it. Specifically, the use of a set of crucial identities, relating the conventional and the 
BFM Green's functions, allow one to convert the new SD series into a dynamical equation 
involving either the conventional or the BFM gluon self-energy only. 

The paper is organized as follows. In Section [III we briefiy review the most salient features 
of the PT methodology and explain qualitatively how the new SD series is obtained and what 
are its main advantages compared to the conventional SD series, focusing on the truncation 
possibilities it offers. In Section IIIII we introduce the notation and the formal machinery 
that will be used in the actual derivation of the PT Schwinger-Dyson equations of QCD. 
We focus particularly on the Batalin-Vilkovisky formalism and the plethora of relations that 
it furnishes for the various fundamental and auxiliary Green's functions appearing in our 
construction. In Section IIVI we present the central result of this work, namely the detailed 
derivation of the new set of SDKs. There are three main subsection, dedicated to the 
construction of the new SDEs for the quark-gluon vertex, the three-gluon vertex, and finally 
the gluon propagator. In Section |V] we discuss some of the main practical implications of 
the new SD series, and present our conclusions and outlook. 

The paper contains five appendices. In Appendix |A] we discuss some subtleties of the 
extension of the PT algorithm beyond one loop, and in particular how to identify unam- 
biguously the subset of three-gluon vertices that must undergo the PT decomposition. In 
Appendix [B] we describe the general strategy for carrying out the renormalization procedure 
to the new SD series obtained within the PT. Finally, the last four appendices furnish the 
derivation of several instrumental formulas employed in Section HVl together with a complete 
set of Feynman rules. 

II. THE NEW SDE SERIES: GENERAL PHILOSOPHY AND MAIN RESULTS 

In this section we present the general ideas and outline the basic philosophy of our 
approach, before diving into the complexities of the full SDE construction. The style of this 
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section is rather qualitative; thus the reader who would like to skip the technicalities can 
get an overview of the theoretical and practical advantages offered by the new SD series, 
compared to the conventional formalism. 



A. The difficulties with the conventional formulation 

Let us first focus on the conventional SD series for the gluon self-energy. Defining the 
transverse projector 

P.M=9.p-^-^, (2.1) 
we have for the full gluon propagator in the Feynman gauge = 1) 



(2.2) 



with A"^(g) = 5°-'^ ^apiq) (in what follows color factors will be omitted whenever possible). 
The scalar function A(g^) is related to the all-order gluon self-energy 

T^aM = Pa/3(g)n(g2), (2.3) 

through 

Since liapiq) has been defined in (12. 4p with the imaginary factor i factored out in front, it 
is simply given by the corresponding Feynman diagrams in Minkowski space. The inverse 
of /S.ap{q) can be found by requiring that 

^^TM{^-n{q) = ^^'9l (2.5) 



and is given by 



or, equivalent ly. 



Kliq) = iPaM^'\q') + n^q^^ (2-6) 



Kliq)=igapq^ -liaM- (2-7) 
In Fig. [1] we show the SDK satisfied by the gluon self-energy. It reads 
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A-i(g2)p^^(g) =g2p^^(g)+Z^(a,)„^. (2.8) 



i=l 



7 



^AA#^AA/ 

Aa Ab 





#\AAA/ 



02 



(«3) 



(04) 




(as) 




FIG. 1: Schwinger-Dyson equation satisfied by the gluon self-energy. The symmetry factors of the 
diagrams are s(ai,a2,«5) = 1/2, 5(03,06) = —1, 5(04) = 1/6. White blobs represent connected 
Green's functions while black blobs represent one particle irreducible ones. 

Of course, in addition to the SDE for the gluon propagator, one must also include the 
corresponding SDE for the vertices; they are normally expressed as skeleton expansions in 
terms of the corresponding connected multi-particle kernels (the proper treatment of the 
vertices is presented in Section [TV]) . 

The main theoretical problem one encounters when dealing with the SDE given above is 
the fact that it cannot be truncated in a physically meaningful way. The most obvious man- 
ifestation of this drawback is the following: after the truncation the fundamental Eq. (11 .ip 
is violated. To recognize the origin of this difficulty, note that Eq. (11. ip translates at the 
level of the SDE to the statement 
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H)af3 



0. 



(2.9) 



i=l 



The diagrammatic verification of (12. 9p . i.e., through contraction of the individual graphs 
by g", is practically very difficult, essentially due to the complicated STIs satisfied by the 



lat of the conventional 



vertices involved. The most typical example of such an STI is t 
three-gluon vertex Va^y{q, ki, ^2) (all momenta entering), given by ,31 1 

q'^T^^^iq, h, h) = [q^D{q)] {A'\kl)PJ {k2)H^^{h, A^a) - A-\kl)P;!{k,)H,^{k2, k^)} , 

(2.10) 

where the auxiliary function Hajs is defined in Fig. [21 The kernel JC appearing in this function 
is the conventional connected ghost-ghost-gluon-gluon kernel appearing in the usual QCD 
skeleton expansion 4, l3^. Notice also that Hai3{k,q) is related to the conventional gluon- 
ghost vertex F^(/c, q) (with k the gluon and q the anti-ghost momentum) by q"'Hai3{k, q) = 



Ha(5{kl, h) = 9al3 + ^ ^ h 



FIG. 2: The auxiliary function H appearing in the three-gluon vertex STL Gray blobs represent 
(connected) Schwinger-Dyson kernels (in this specific case the ghost-gluon kernel K- appearing in 
the usual QCD skeleton expansion). 



Tp{k,q) y,l3l|,l32|,l33|. 



In addition, some of the pertinent STIs are either too complicated, such as that of the 
conventional four-gluon vertex, or they cannot be cast in a particularly convenient form. 
For instance, in the case of the conventional gluon-ghost vertex, F^(g,p), the STI that one 
may obtain formally for q^T^{q,p) is the sum of two terms one of which is p^T ^[q,p)\ this 
clearly limits its usefulness in applications. 

The main practical consequence of these complicated STIs is that one cannot truncate 
(12. 8p in any obvious way without violating the transversality of the resulting liapiq)- For 
example, keeping only graphs (ai) and (02) is not correct even perturbatively, since the 
ghost loop is crucial for the transversality of IIq,^ already at one-loop; adding (03) is still not 
sufficient for a SD analysis, because (beyond one-loop) g"[(ai) + (02) + {a^)]af3 7^ 0. 



B. The pinch technique 

The FT jl^ . Il^ is a particular algorithm for rearranging the perturbative series in such a 
way as to obtain new Green's functions that are independent of the gauge-fixing parameter, 
and satisfy to all orders ghost-free, QED-like Wis, instead of the usual STIs. The original 
motivation for resorting to it was precisely to devise a truncation scheme for the SDE that 
would preserve manifestly the gauge-invariance of the answer at every step. 

Let us emphasize from the beginning that, to date, there is no formal definition of the 
FT procedure at the level of the functional integral defining the theory. In particular, let 
us assume that the path integral has been defined using an arbitrary gauge-fixing procedure 
{i.e. linear covariant gauges); then, there is no known a priori procedure (such as, e.g., 
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functional differentiation with respect to some combination of appropriately defined sources) 
that would furnish directly the gauge- independent PT Green's functions. The definition of 
the PT procedure is operational, and is intimately linked to the diagrammatic expansion of 
the theory {i.e. one must know the Feynman rules). In fact, the starting point of the PT 
construction can be any gauge-fixing scheme that furnishes a set of well-defined Feynman 
rules and gauge-independent physical observables. Specifically, one operates at a certain 
well-defined subset of diagrams, and the subsequent rearrangements give rise to the same 
gfp-independent PT answer, regardless of the gauge-fixing scheme chosen for deriving the 
Feynman rules. Note however that, as the present paper amply demonstrates (and as has 
already been emphasized in some of the earlier literature), the PT is not diagrammatic, in 
the sense that one does not need to operate on individual graphs but rather on a handful of 
classes of diagrams (each one containing an infinite number of individual graphs). This is 
the enormous advantage of formulating the PT at the SD level. 

The aforementioned rearrangements of the PT are collectively implemented through the 
systematic use of the STIs satisfied by certain Green's functions and kernels; the latter 
constitute standard ingredients in the ordinary perturbative expansion or the SDK of the 
various n-point Green's functions. In the Feynman gauge, which is by far the most convenient 
choice, the relevant STIs are triggered by the action of a very special set of longitudinal 
momenta. Specifically, consider the subset of Feynman diagrams that have at least one 
external three-gluon vertex 



By "external" we mean a vertex that has one of its legs irrigated by a physical momen- 
tum, to be denoted by q, as opposed to a virtual momentum, that is being integrated over 
in the Feynman graph [30] (a detailed discussion of why only external vertices can pinch 
while all other three-gluon vertices inside the loops should remain unchanged is provided in 
Appendix [All . Then Taf^uiq, ki, is decomposed as jl^ 



Evidently the above decomposition assigns a special role to the leg carrying the physical 




(2.11) 






(2.12) 



10 



momentum q, and allows r^^^(g, ki, /C2) to satisfy the tree- level WI 

g-r^^,(g, h, h) = {kl - kl)g,,, (2.13) 

where the rhs is the difference of two (tree-level) inverse propagators in the Feynman gauge. 
Note that T^^^^iq, ki, /C2) (i) is Bose-symmetric only with respect to the fi and z/ legs, and (ii) 
it coincides with the BFM three-gluon vertex involving a background gluon, Aa{q), and two 
quantum gluons, A^{ki) and Ay{k2), in the Feynman gauge {i.e., when the quantum gauge- 
fixing parameter, ^q, is chosen to be = 1). The term r^^^,(g, fci, k2), which in configuration 
space corresponds to a pure divergence, plays the central role in the PT construction; indeed, 
the main thrust of most PT demonstrations (in this article and many others before) is to 
essentially track down the precise action of the momenta contained inside P^^^. Specifically, 
^o^i/ contains the longitudinal "pinching" momenta, which will get contracted with the 
kernels and Green's functions nested inside the remaining part of the diagram, triggering 
the corresponding STIs; this, in turn, will produce the highly non-trivial rearrangements 
of the various terms characteristic of the PT. Quite remarkably, all these rearrangements 
finally amount to the modification of the ghost sector of the theory, reproducing dynamically 
the corresponding ghost sector of the BFM, leaving no residual terms behind. 

The simplest example that demonstrates the action of the pinching momenta is the one- 
loop construction of the PT gluon self-energy (see subsection III D^ : the STI triggered inside 
the conventional one- loop diagram (a) in Figl3]is simply the tree- level version of Eq. fl2.10l) . 
namely 

fcfr„^.(g, k^) = q^Pauiq) - klP^uih), (2.14) 
fc2"r<,^.(g, K k2) = klPo^^ih) - q^PaM)- (2-15) 

The terms proportional to an inverse propagator of the external leg {i.e., to g^), will cancel 
(when embedded into a physical process!) against similar contributions from other graphs 
{e.g., vertex-graphs), also produced by the corresponding action of the pinching momenta 
inside them; in this case the pinching momenta literary "pinch out" internal quark lines 
(hence the name of the technique). Thus, effectively, inside a physical process, this particular 
subset of pinching contributions can be discarded altogether; this is the shortcut introduced 
in the "intrinsic" |l^. The remaining pinching terms, namely those proportional to the 
internal leg, are instrumental for obtaining the PT answer; in particular, they symmetrize 
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the original ghost sector, so that it will finally coincide with the BFM ghost sector at 
= 1 (see again subsection III D|) . 

C. Pinch technique and background field method: some conceptual issues 

As mentioned in the Introduction, the (gauge-independent) PT n-point functions coincide 
with the BFM n-point functions when the latter are computed at the special value = 1 
(BFG). Even though this correspondence (and its correct interpretation) has been addressed 
in various places in the literature, it may be useful to present a brief overview of some of 
the main subtleties associated with it. 

i. The objective of the PT construction is not to derive diagrammatically the BFG, but 
rather to exploit the underlying BRST symmetry in order to expose a large number 
of cancellations, and eventually define gauge-independent Green's functions satisfying 
abelian Wis. In fact, it was after more than a decade of independent PT activity 
(when practically all one loop calculations had been carried out both in QCD and the 
electroweak sector) when the aforementioned correspondence was discovered (i.e. the 
PT results already existed, and then it was realized that they coincide with the results 
of the BFG). Thus, while it is a remarkable and extraordinarily useful result that the 
PT Green's functions can also be calculated in the BFG, this needs a very extensive 
demonstration. Therefore, the correspondence must be verified at the end of the PT 
construction and should not be assumed beforehand. 

a. It is well known that, at any order, the S'-matrix S of the BFM, is equal to that in the 
conventional linear (R^) gauges, i.e., S = S There is no way, however, to deduce from 
this equality the PT-BFG correspondence. Writing S = TAT + B and S = TAT + B, 
using that the box diagrams are equal in both schemes, i.e., B = B, and, finally, 
observing that the PT does not change the unique S'-matrix, one can deduce that 
PAP = PAP, and hence that PAP = PAP. But from this does not follow that A = A 
nor that P = P; one must prove explicitly the equality for individual Green's functions. 

Hi. We emphasize that the PT is a way of enforcing gauge independence (and several other 
physical properties) for off-shell Green's functions; the BFM, in a general gauge, is 
not. This is refiected in the fact that the BFM n-point functions are gauge- invariant, 
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in the sense that they satisfy (by construction) QED-hke Wis, but are not gauge- 
independent, i.e., they depend exphcitly on ^q. Had the BFM n-point functions 
been ^g-independent, in addition to being gauge- invariant, there would be no need for 
introducing independently the PT. 

iv. Notice that the ^g-dependent BFM Green's functions are not physically equivalent. 
This is best seen in theories with spontaneous symmetry breaking: the dependence 
of the BFM Green's functions on C,q gives rise to unphysical thresholds inside these 
Green's functions for 7^ 1, a fact which limits their usefulness for resummation 



purposes 



2l| . Only the case of the BFG is free from unphysical poles; that's because 



then (and only then) the BFM results collapse to the physical PT Green's functions. 

V. The PT procedure has no a-priori knowledge of the BFM built into it, despite the 
fact that the splitting of the regular three gluon vertex given in (I2.12p suggests such 
a preference. Indeed, while coincides with the BFG vertex, it only furnishes one 
piece of the final answer. As already mentioned, the non-trivial part of the PT con- 
struction resides in what happens when the F^ part of the vertex gets contracted with 
the Green's functions and kernels nested inside the corresponding SD diagram. The 
correspondence with the BFG works finally only because the WI triggered by F^ con- 
spire in such a way as to reproduce dynamically the BFM ghost sector at C,q = 1, and 
nothing more. There is no a-priori way of knowing that this will indeed happen, and 
hence the need for the detailed demonstration presented in the next sections. 

vi. Amplifying the previous point, notice that the PT works perfectly well in the context 
of non-covariant gauges (in fact it was first carried out in such a gauge), where the 
ghosts are decoupled from the S'-matrix. Spectacularly enough, the PT procedure 
produces completely dynamically the necessary ghost sector, from the STIs that are 
triggered. 

vii. Perhaps the most compelling fact that demonstrates that the PT and the BFM are 
intrinsically two completely disparate methods is the following: one can apply the PT 
within the BFM. For example, the PT can be used to combine pieces of Feynman 
graphs in the background Landau gauge, just as in any other gauge, and the usual PT 
results (those of the BFG) emerge. Operationally this is easy to understand: away 
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from = 1 even in the BFM there are longitudinal (pinching momenta) that will 
initiate the pinching procedure. Ultimately, the BFG is singled out because of the 
total absence, in this particular gauge, of any such longitudinal momenta. 

via. We emphasize that the PT construction goes through unaltered under circumstances 
where the BFM Feynman rules cannot even be applied. Specifically, if instead of 
an S"- matrix element one were to consider a different observable, such as a current- 
current correlation function or a Wilson loop (as was in fact done by Cornwall in the 
original formulation [l^, and more recently in 2^) one could not start out using the 
background Feynman rules, because all fields appearing inside the first non-trivial loop 
are quantum ones. Instead, by following the PT rearrangement inside these physical 
amplitudes the unique PT answer emerges again. 



D. The pinch technique as a gauge-invariant truncation scheme 

Let us now see how the standard one-loop PT construction contains the seed of a gauge- 
invariant truncation scheme for the SDK of the gluon self-energy. This exercise may seem 
trivial at first, in the sense that no truncation is really needed, given that the two diagrams 
comprising the full answer are elementary to calculate. However, it illustrates exactly how 
the PT rearrangement furnishes a transverse one-loop approximation for the conventional 
gluon self-energy, even if the (modified) ghost loop is omitted. 

In what follows we will use dimensional regularization, and will employ the short-hand 
notation J^, = J d'^k/{2TTY, where c? = 4 — e is the space-time dimension. The conventional 
one-loop self-energy in the Feynman gauge, to be denoted by Il^^p{q), is given by the diagrams 
(a) and (6) in Fig. [3] (we set the "seaguH"-type contributions directly to zero, using the 
standard result Jj^k'"^ = 0). As is well known, neither (a) nor (6) is transverse, and it is 
only their sum that furnishes a transverse answer for Il^^^{q). Specifically, setting 

with Ca the Casimir eigenvalue of the adjoint representation {Ca = N for SU{N)), and 
dropping irrelevant constants, we have 

{a)a/3 = ^/(g^)(19gV/3 - 22g„g^), 
14 



PT 



A a 




4* 



FIG. 3: The conventional one-loop gluon self-energy before (first line) and after (second line) the 
PT rearrangement. A gray circle at the end of an external gluon line denotes that the corresponding 
gluon behaves as if it were a background gluon. 



r(i) 



f{(l'^){(1^9ai3 + 2qaqi3), 



Kil) = 5q'f{q')PM- 



(2.17) 



The application of the PT amounts to carrying out the following rearrangement of the 
two elementary three-gluon vertices 



afiu] 



^ (3 +^/3 



(2.18) 



Then, using the elementary Wis of Eq.s fl2.14|) and (12.151) we have that 



pP pP 



2kakp + {k^qp + qakp), 



(2.19) 
(2.20) 



where some terms have been set to zero by virtue of the dimensional regularization result 
A;~^ = 0. Thus, one can cast n^^j(g) in the following form: 



pF pF/i!/ 
k^{k + qy 



{2k + q)a{2k + q)^ 
k'^ik + g)2 



- 2CV 



k^{k + qy 



(2.21) 



It is elementary to verify that each of the two terms in the square bracket on the rhs of (I2.2ip 
are transverse; thus the PT rearrangement has created three manifestly transverse structures. 
That in itself might not be so important, if it were not for the fact that these structures 
admit a special diagrammatic representation and a unique field-theoretic interpretation. 
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FIG. 4: Schematic representation of the pinching contributions one needs to consider away from 
the Feynman gauge ^ = 1. 



Specifically, the two terms in the square bracket correspond precisely to diagrams (a) and 
(6) defining the one-loop gluon self-energy in the BFG, to be denoted by n|^^j(g)|^Q=i; note 
in particular the symmetrized gluon-ghost coupling. Thus, it is as if the external gluons 
in (a) and (b) had been converted dynamically into background ones. As explained in jl3| . 
the third term on the rhs of (12.211) is the one-loop expression of a special auxiliary Green's 
function, to be defined shortly; it corresponds to diagram (c) in Fig. [3], and is generated 
from the first term on the rhs of fl2.14p . The one- loop PT self-energy, to be denoted by 
n^^j(g), is obtaine d by simply dropping this last term from the rhs of fl2.2ip : this defines 
the "intrinsic" PT ISj. The completely equivalent way of saying this, corresponding to the 
"S'-matrix" PT [l4], is that the term corresponding to graph (c) cancels exactly against 
a propagator-like contribution extracted from the vertex graphs contributing to the full S'- 
matrix element that one considers. Notice that this is true only in the Feynman gauge: away 
from ^ = 1 additional pinching contributions need to be considered, e.g., the ones coming 
from box and self-energy correction diagrams (Fig. Hj). 

To see how the aforementioned cancellation comes about, let us now imagine that the 
gluon self-energy 11^^^ (g) is embedded into a physical process, such as the S'-matrix element 
for the quark-quark elastic scattering process ip{ri)ilj{r2) —>■ ^(pi)V'(P2); with q = ri — r2 = 
P2 — Pi being the momentum transfer. The one-loop quark-gluon vertex consists of the two 
graphs shown in Fig. [5l Let us concentrate on the non-Abelian diagram [d), and carry out 
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^{Pl) 



FIG. 5: The one-loop quark-gluon vertex appearing in the quark-quark elastic scattering process. 



the vertex decomposition of Eq. (12.121) [l5l | : 



{d) 



^ -Ik 



k\k + qy 



k^{k + qY 



(2.22) 



The first term on the rhs of the second line is the pure vertex-like part of (rf), while the 
second term is purely propagator-like, as can be easily established using the elementary WI 



Kl" = {f^ + i) - m) - - m) 



(2.23) 



The first term on the rhs of (I2.23P removes (pinches out) the internal bare quark propagator 
'S'o(^ + i>), whereas the second vanishes on shell, since u{p2){,'^2 — 'm) = and (^i — m)u(pi) = 
0. Thus, 



2i 



k^ik + qy - J^k\k + qf^''- 

The self-energy-like contribution from the two vertex graphs (mirror graph included), to be 
denoted by n^^j^(g), is given by (longitudinal pieces may be added for free, due to current 
conservation) 



(2.24) 



ni7(g) = 2CV/ 

Jk 



k^{k + qy' 

The PT one-loop quark-gluon vertex, to be denoted by r^(pi,p2), is given by [l^ 



(2.25) 



2,a 



k k^k + qf 

rs^''KH^2)rs(-^\lk+^,)j^ 
k^ 



(2.26) 
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where C/ is the Casimir eigenvalue of the fundamental representation [C/ = (A^^ — 1)/2N 
for SU{N)]. Now it is easy to derive the QED-like WI that the T'^{pi,P2) satisfies. Using 
(I2.13p . we have that 

Jk J 

(2.27) 



where S(^) is the one-loop quark self-energy in the Feynman gauge [20|, |2 
Returning to the gluon self-energy, n|^^^(g) is defined as 

nS(^) = n2(?) + nS^)(g). 



(2.28) 



After carrying out the integrals one obtains 



and thus [1 



(c)a/3 



-(^q'fiq')Pa^{q), 



= nq'f{q')PM=^%)k,=i- 



(2.29) 



(2.30) 



Note that the second line of fl2.30p expresses the PT-BFG correspondence at one loop [26 1. 
Then, Eq. f l2.2ip assumes the alternative form 



Kiiq) = Kkq) + ic 



Jap- 



(2.31) 



As has been explained in detail in 



34| . and as we will see in the following sections, Eq. fl2.3ip 



is the one-loop version of a general identity 35|, which we will call "Background-Quantum" 
identity (BQl) [34j], relating the conventional and the BFM self-energies in terms of an 
auxiliary Green's function, corresponding to graph (c). This identity is valid to all orders in 
perturbation theory, as well as non-perturbatively, and may be obtained either formally, by 
resorting to the Batalin-Vilkovisky (BV) formalism, or diagrammatically, as a by-product 
of the FT rearrangement of the conventional SD series; as we will see, in this latter case no 
reference to the BV formalism is necessary. 
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FIG. 6: The new Schwinger-Dyson series projected out dynamically by the PT algorithm. The 
symmetry factors are in this case 5(^1,^2,^6) = 1/2, 5(^5) = 1/6, and all the remaining diagrams 
have s = —1. 

Let us now focus on n^^^(g) and imagine for a moment that no ghost loops may be 
considered when computing it, i.e., the graphs must be omitted; in a SDE context this 
"omission" would amount to a "truncation" of the series. One may still obtain a transverse 
approximation for 11^^^ (g) with no ghost-loop, given by 

nS(g) = {a)ap + {cU = Wf{q^)PaM- (2.32) 

Interestingly enough, the PT rearrangement offers already at one-loop the ability to truncate 
gauge-invariantly, i.e., preserving the transversality of the truncated answer. 



E. The new Schwinger-Dyson series 

The implementation of the PT at the level of the SDE has been studied first in the 

n 

context of scalar QED [36]. The corresponding construction in the case of quarkless QCD 
has been recently carried out in a short communication [13], where we restricted ourselves 
to the SDE of the gluon propagator. As has been explained there, the PT rearrangement 
gives rise dynamically to a new SD series (see Fig. E]), with the following characteristics: 
on the rhs we have graphs that are made out of new vertices, but contain inside them the 
same gluon propagator as before, namely A„/3(g). The new vertices, to be denoted by P"^", 
-panm^ F^™"^, P^™"''', correspoud precisely to the Feynman rules of the BFM in the Feynman 
gauge, i.e., as already seen explicitly in the one-loop case, it is as if the external gluon had 
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FIG. 7: Diagrammatic representation of the auxiliary functions H and A. 

been converted dynamically into a background gluon. The Ihs is composed from the sum of 
three terms: in addition to the term A~^{q'^)Pai3{q), present there from the beginning, we 
have two additional contributions, 2G{q'^)A~^{q'^)Pai3{q) and G"^ [q"^) A'^ {q'^) Papiq) , which 
appear during the PT rearrangement of the rhs (and are subsequently carried to the Ihs). 
The quantity G{q'^) is a special function, defined in terms of the gluon and ghost propagators 
as well as the auxiliary function of Fig. [2l Specifically, define the following two-point 
function Aa/siq), (we suppress color indices) 

^afsiq) = Ca [ Hf^D{k)A>^^{q - k) H,p{q - k, -q), (2.33) 
Jk 

with the diagrammatic representation shown in Fig. [71 Then, G{q'^) is defined as i times the 
component of lyap^q) multiplying Qap, namely 

KM = igc.pG{q^) + . . . , (2.34) 

where the omitted terms are proportional to qaq/i- Thus, the term appearing on the Ihs of 
the new SDE is A^-'^(g^)[l + G{q^)Y Paii{q) . So, one may write schematically 

11 

A-\q')[l + G{q^)fP^p{q) = q'P^ + ^J2^d^U^ (2-35) 

1=1 

or, equivalently, casting it into a more conventional form with the inverse of the unknown 
quantity isolated on the Ihs, as 

A-'tf )P„M<!) = ''"'"'ilV^Gpl'''''''"' - (2-36) 

This new SD series has a very special structure. Let us first separate the diagrams on 
the rhs into four obvious categories: one-loop (dressed) gluonic contributions [{di) and (^2)], 
one-loop ghost contributions [(^3) and (d^)], two-loop gluonic contributions [(^5) and (c/g)], 
and two-loop ghost contributions [{dj), (rfg), {dg) and (rfio)]. It turns out that, by virtue 
of the all-order WI satisfied by the full vertices F^^", F^""", T'^J^^, ri^^""" appearing in the 
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various diagrams, the contribution of each of the four subgroups is individually transverse. 
Specifically, the four fundamental all-order Wis are given by 

?"r— (g, h, h) = [A-l{k,) - A-Jih)] , (2.37) 

g"fr'"(g, h, h) = ign^ [D-\k,) - D-\h)\ , (2.38) 
?°r~(g, fci, h) = gr''^V%^{q + h) + gr'^VtZiq + h) 

+ ^7r'^"'rJ;(g + A;i,fc2,A;3), (2.39) 
g"f~(g, k,, A;2, k^) = -gf-^^V^-^iq + k,, k^, k,) - (?/'^"T--(A;i, q + k^, k,) 

- gr^r^^'ik^k^q + k,). (2.40) 



Using these Wis one may show after some elementary operations that 



37| 



[(^i) + {d2)U = 0' 

lids) + {d,)U = 0, 
[(4) + {de)U = 0' 
q" [{dr) + ids) + (dg) + {dio)U = ^- (2-41) 

[Notice that the one-loop dressed fermionic contributions (dn) trivially satisfy this transver- 
sality property.] 

As has been pointed out in [l3^, this special property has far-reaching practical conse- 
quences for the treatment of the SD series. Specifically, it furnishes a systematic truncation 
scheme that preserves the transversality of the answer. For example, keeping only the dia- 
grams in the first group, we obtain the truncated SDE 

A-lr^2^p _ q''Pap{q) + i[{di) + {d2)U 

and from the first equation of fl2.4ip we know that [{di) + {d2)]ai3 is transverse, i.e., 
[{di) + {d2)]ai3 = {d — l)^'^[{di) + {d2)]'^Paf3{q) ■ Thus, the transverse projector Papiq) ap- 
pears exactly on both sides of fl2.42p : one may subsequently isolate the scalar cofactors on 
both sides obtaining a scalar equation of the form 

[q) [1 + G(g2)]2 ■ ^^■'''^) 

A truncated equation similar to fl2.42p may be written for any other of the four groups, or 
for sums of these groups, without compromising the transversality of the answer. The price 
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one has to pay for this advantageous situation is that one must consider in addition the 
equation determining G'(g^), i.e., the Qa/s part of Eq. fl2.33p . This price is, however, rather 
modest, given that Eq. fl2.33p may be approximated introducing, for example, a dressed- 
loop expansion (see Fig 12]), without jeopardizing the transversahty of Uajsiq), given that 
[1 + G'(g^)]^ affects only the size of the scalar prefactor. 

In going from Eq. fl2.35p to Eq. (12.360 one essentially chooses to retain the original 
propagator A{q) as the unknown quantity, to be dynamically determined from the SDE. 
There is, of course, an alternative strategy: one may define a new "variable" from the 
quantity appearing on the Ihs (12.351) . namely 

A{q)^[l + G{q')]-'A{q), (2.44) 

which leads to a new form for (I2.35p . 

11 

A-\q^)P^p{q) = q^Pap{q)+^Y.^d,)^p. (2.45) 

i=l 

Obviously, the special transversahty properties established above holds as well for Eq. (12.451) ; 
for example, one may truncate it gauge-invariantly as 

A-\q^)P^^{q) = q^P^M + i[{d{) + {d^)^. (2.46) 

Should one opt for treating A(g) as the new unknown quantity, then an additional step 
must be carried out: one must use (12.441) to rewrite the entire rhs of (I2.45P in terms of A 
instead of A, i.e., carry out the replacement A — [1 + A inside every diagram on the 
rhs of Eq. (I2.45P that contains A's. 

Let us discuss further these two versions of the SDE. Eq. (I2.42p furnishes a gauge-invariant 
approximation for the conventional gluon self-energy A(g), whereas Eq. (12.450 is the gauge- 
invariant approximation for the effective PT self-energy A. The crucial point is that one 
may switch from one to the other by means of Eq. (12.440 . For practical purposes this means 
for example, that one may get a gauge-invariant approximation not just for the PT quantity 
(background Feynman gauge) but also for the conventional self-energy computed in the 
Feynman gauge. Eq. ( 12.440 . which is the all-order generalization of the one-loop relation 
given in Eq. (12.311) . plays an instrumental role in this entire construction, allowing one to 
convert the SDE series into a dynamical equation for either A(g) or A(g). 
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III. THE FORMAL MACHINERY 



The extension of the PT algorithm to the SDEs of QCD is a challenging exercise, mainly 
due to the large amount of different Green's functions one needs to manipulate in the process. 
Most of these Green's functions are generated when longitudinal momenta trigger the STIs 
satisfied by specific subsets of fully dressed vertices appearing in the ordinary perturbative 
expansion. Due to the non-linearity of the BRST transformation [see Eq. (13.81) below] they 
involve composite operators; specifically, they are of the type {0\T[s^{x) ■ ■ -JlO) with s the 
BRST operator and $ a generic QCD field. 

It turns out that the most efficient framework for dealing with this type of quantities is 
the BV formalism, allowing the construction of these auxiliary (ghost) Green's functions in 
terms of a well-defined set of Feynman rules. In addition, this formalism furnishes a set of 
useful identities (the BQIs mentioned in the previous section), relating Green's functions 
involving background fields to Green's functions involving quantum fields. 

In this section, we fix our conventions, present the QCD Lagrangian and its gauge- fixing 
procedure (concentrating, in particular, on the conventional gauges and the BFM), and 
briefly review the BV formalism. Then, we proceed to describe how one can extract from 
the master equations the all-order STIs and BQIs needed in the coming PT construction, 
postponing their actual derivation to the Appendix [D] and [El We will also describe how to 
derive the so-called Faddeev- Popov equations (see also Appendix [C]) . Finally, in the process 
of describing all the above topics, we will introduce a particularly compact notation for 
Green's functions, which encodes unambiguously all relevant information {i.e., the particle 
content, Lorentz and color structure, and momenta flow). 

A. QCD Lagrangian and gauge fixing schemes 

Throughout the paper we will adopt the conventions of the book by Peskin & Schroder 



381]. The QCD Lagrangian density is given by 

C = Ci + Cgf + Cfpg- (3.1) 
£i represents the gauge invariant SU (3) Lagrangian, namely 

£i = -^Fr^;, + {zYV, - m),^. ^1, (3.2) 
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where a — 1, . . . , 8 (respectively i,j — 1, 2, 3) is the color index for the adjoint (respectively 
fundamental) representation, while "f" represents the flavor index. The field strength is 

F;, = d,At - d^Al + gr'^A^Al, (3.3) 

and the covariant derivative is defined according to 

{V,),^ = d,{l),^ - igA^n^, (3.4) 

with g the (strong) coupling constant. Finally, the SU{3) generators t"' satisfy the commu- 
tation relations 

[r, t^] = ir^H", (3.5) 
with f"'^'^ the totally antisymmetric SU (3) structure functions. 

£gf and >Cfpg represent respectively the (covariant) gauge fixing Lagrangian and its asso- 
ciated Faddeev-Popov ghost term. The most general way of writing these terms is through 
the expressions 

Upg = -c«s^«. (3.7) 

In the formulas above is the gauge fixing function, and the are auxiliary, non- 
dynamical fields (the so called Nakanishi-Lautrup multipliers) that can be eliminated 
through their (trivial) equations of motion; c" (respectively, c") are the ghost (respectively, 
anti-ghost) fields, and, finally, s is the BRST operator, with the BRST transformations of 
the QCD fields given by 

sAl = a^c« + gr^'Alc'^ sc- = -^^/"''=c^c^ 
stlj\ = igc''{r)ijilji se = B% 

= -igcyi{t'')ji sB" = 0. (3.8) 

We thus see that the sum of the gauge fixing and Faddev- Popov terms can be written as a 
total BRST variation: 

Cgf + Cfpg = s (^c'^.F" - |c"S"^ . (3.9) 

This is of course expected, since it is well known that total BRST variations cannot appear in 
the physical spectrum of the theory. For our purposes, the gauge-fixing functions of interest 
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are the ones corresponding to the (renormahzable ^ gauges) and the BFM, which we 
describe in what follows. 



1. In the usual gauges, the gauge fixing function is chosen to be J^^^ — d'^A'l^; therefore 
one finds 

Cgf = (3.10) 

£fpg = a^c'^a^c" + ^r^<^(a'^c")4c^ (3.11) 

2. In the case of the BFM, one starts by splitting the gluon field into a background part, 
A^, and a quantum part, A^. Notice that the ERST variation of the background field 
will be zero, but the latter will enter in the variation of the quantum one, since 

sAl = d,c'^ + gr'^{A'^ + A'^y. (3.12) 

The gauge fixing function is 

BFM — \^ ^t^J 

cabc 'Ab 



d^Al + gr'^A^A'::, (3.13) 



which gives in turn 



£gf = Yi{d>^Alf + -^gr''^{d>^Al)AlA: + ^^g' P'^f^^AlAlAlA^,, (3.14) 
£ppG = di'^d^c" + gf''^''{d>'e)Alc^ + gr^%d''^)A^^c^ - g^'^'c'^A^^id^'c^) 

- gT'^r'^c'^A'^iA^^ + Apc'^. (3.15) 

We thus see the appearance of the characteristic ghost sector for the interaction with 
background gluons, consisting in a symmetric Acc ghost vertex and a four particle 
A Ace one. 

B. Green's functions: conventions 

The Green's functions of the theory can be constructed in terms of timc-ordcrcd products 
of free fields $5 ' ' ' vertices of the interaction Lagrangian >Cint (constructed from the 
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pieces of L which are not bihnear in the fields) through the standard Gell-Man-Low formula 
for the IPI truncated Green's functions 

r$,...$„(xi,...,x„) = (T[$i(xi)---$„(x„)])iP^ 

= (T[$?(xi)---$°(xO]exp(-z y"d^xA„0)^P^ (3.16) 

The complete set of Green's functions can be handled most efficiently by introducing a 
generating functional, which in Fourier space reads 

rt*^] = E ^ / n ^'^^ ^^^V,)^l{Vx) ■ ■ ■ $n(Pn)r*,...*„(Pl, . . . ,Pn), (3.17) 
n=0 ■ i=0 j=l 

with Pi the (in-going) momentum of the field. Since in perturbation theory r<j,^...<j,^ is a 
formal power series in h, we will denote its m-loop contribution as rL .^..(j, . Then, in terms 
of the generating functional r[$] any of the Green's function of the theory can be obtained 
by means of functional derivatives: 



(3.18) 

*i=0 



where $(p) denotes the Fourier transform of $(x) and our convention on the external mo- 
menta is summarized in Fig. [HI From the definition given in Eq. f l3.18p it follows that the 
Green's functions z~"r,j,j...$^ are simply given by the corresponding Feynman diagrams in 
Minkowski space. Finally, notice that upon inversion of two (adjacent) fields we have 

r$i...$,$,+i...$„(pi, . . . ,Pi,Pi+l, . . . ,Pn) = ±r<i,^...$^_^^$^...$„(pi, . . . ,Pi+l,Pi, . . . ,Pn), (3.19) 

with the minus appearing only when both fields $i and obey Fermi statistics. 

The Green's functions constructed so far are sufficient for building all possible amplitudes 
involved in the S'-matrix computation; however, due to the non-linearity of the BRST trans- 
formations [Eq. (13.81) ]. they do not cover the complete set of Green's functions appearing 
in the STls of the theory (and therefore needed for its renormalization, as well as the PT 
construction). 



C. A brief introduction to the Batalin-Vilkovisky formalism 



In this subsection we review briefiy the BV formalism 39|], which allows one to obtain 
both the STIs as well as the BQIs of the theory at hand. In order to simplify the notation 
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FIG. 8: Our conventions for the (IPI) Green's functions T^-i^...^^{pi, . . . ,pn). All momenta 
P2t ■ ■ ,Pn are assumed to be incoming, and are assigned to the corresponding fields starting from 
the rightmost one. The momentum of the leftmost field <I>i is determined through momentum 
conservation {J2iPi = 0) and will be suppressed. 
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TABLE I: Ghost charges and statistics (B for Bose, F for Fermi) of the QCD fields, anti-fields and 
BFM sources. 

(and since they will not play any role in what follows) we will suppress from now on all 
spinor indices (both flavor and color). 

Let us then start by introducing for each field $ appearing in the theory a correspond- 
ing anti-field, to be denoted by $*. The anti- field $* has opposite statistics with respect 
to $; its ghost charge, gh($*), is related to the ghost charge gh($) of the field $ by 
gh($*) = — 1 — gh($). For convenience, we summarize the ghost charges and statistics of 
the various QCD fields and anti-fields in Tabled Next, we add to the original gauge invari- 
ant Lagrangian a term coupling the anti-fields with the BRST variation of the corresponding 
fields, to get 

^BY = A + -^BRST, 
J-'BRST = 

= /i;'^(9^c" + gT^^A^c^) - -gr^^d^'c^c" + ig^^d^f^jj - igc^'tpetp* . (3.20) 
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Then, the action r'^'')[$, $*] constructed from Cbv, will satisfy the master equation 



0. (3.21) 



To verify this, observe that, on one hand, the terms in bV^^^ jb^ that are independent from 
the anti-fields $* are zero due the BRST (actually the gauge) invariance of the action 

on the other hand, terms in bV'^^^ jb^ that are linear in the anti-fields vanish due to the 
nihilpotency of the BRST operator 

Now, since the anti-fields are external sources, we must constrain them to suitable values 
before we can use the action Y^^^ in calculations of ^-matrix elements 25j. To that end, we 
introduce an arbitrary fermionic functional, with ghost charge -1, and set for all the 

anti-fields 

$* = (3.24) 

Then the action becomes 

r(o)[$,(5^/5$] = r{°^[$] + (5$) 



5$ 

= ri°^[<l>] + s^[<l>], (3.25) 
and therefore, choosing the functional to satisfy the relation 

S^ = jd^x{CGF + CFPG), (3.26) 

we see that the action F'^^^ (obtained from £bv) is equivalent to the gauge-fixed action 
obtained from the original Lagrangian C of Eq. (13.11) . The functional ^' is often referred to 
as the "gauge fixing fermion" . 

The BRST symmetry is crucial for endowing a theory with a unitary S'-matrix and gauge- 
independent physical observables; therefore, it must be implemented to all orders. For doing 
so we establish the quantum corrected version of the master equation (]3.2ip in the form of 
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the STI functional 

6T 6r 



Sim] 



5$* 5$ 

, ST ST bT ST ST ST ST ST ST 



= 0, (3.27) 

where r[$, $*] is now the effective action. 

In order to simphfy the structure of the STI generating functional of Eq. (I3.27p . let 
us notice that the anti-ghost c" and the multiplier have linear BRST transformations; 
therefore they do not present the usual complications (due to non-linearity) of the other QCD 
fields. Together with their corresponding anti-field, they enter bi-linearly in the action, and 
one can write the complete action (which we now explicitly indicate it with a C subscript) 
as a sum of a minimal and non-minimal sector 

rj?^ [$, $*] = r(°) [A"^, ^\ ^, ^\ c^ c*"] + e'^B". (3.28) 

The last term has no effect on the master equation (13.211) . which is satisfied by T^'^^ alone; 
the fields {^^, A*^, ip, ip*, ip, ip*, c", c*"} are then often called minimal variables while and 
5" are referred to as non-minimal variables or "trivial pairs". Equivalently one can intro- 
duce the minimal (or reduced) action by subtracting from the complete one the local term 
corresponding to the gauge-fixing Lagrangian, i.e., 



T = Tc- J d'^xCcF- (3.29) 

In either cases, the result is that the STI functional is now written as 

o/T.xr^n f ,4 f '^r 6T 6T 6T 6T 6T 6T 6T ] , , 

^ ^ J \5A*r^5A'^ 5c*"'Sc"' Sip*6ip Sip 6ip* ) 

In practice, the STIs generated from the functional of Eq. (13.301) coincide with the one 
obtained by the complete one after the implementation of the Faddeev-Popov equation de- 
scribed in the next subsection [4^. One should also keep in mind that the Green's functions 
involving unphysical fields generated by the minimal functional coincide with the ones gener- 
ated by the complete functional only up to constant terms proportional to the gauge fixing 
parameter, e.g., Tai,a^{(i) = T^ A^il) " '^^^^Qfj.^u- We will discuss further the differences 
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between employing the complete and minimal generating functionals in the Appendix [D] 
andm 

Taking functional derivatives of iS(r)[$] and setting afterwards all fields and anti- fields 
to zero will generate the complete set of the all-order STIs of the theory; this is in exact 
analogy to what happens with the effective action, where taking functional derivatives of 
r[$] and setting afterwards all fields to zero generates the Green's functions of the theory, 
see Eq. (13.181) . However, in order to reach meaningful expressions, one needs to keep in mind 
that: 

1. S(T) has ghost charge 1; 

2. functions with non-zero ghost charge vanish, since the ghost charge is a conserved 
quantity. 

Thus, in order to extract non-zero identities from Eq. (13.301) one needs to differentiate the 
latter with respect to a combination of fields, containing either one ghost field, or two ghost 
fields and one anti-field. The only exception to this rule is when differentiating with respect 
to a ghost anti-field, which needs to be compensated by three ghost fields. In particular, 
identities involving one or more gauge fields are obtained by differentiating Eq. (I3.30p with 
respect to the set of fields in which one gauge boson has been replaced by the corresponding 
ghost field. This is due to the fact that the linear part of the BRST transformation of the 
gauge field is proportional to the ghost field: sA^lunear = <9^c". For completeness we notice 
that, for obtaining STIs involving Green's functions that contain ghost fields, one ghost 
field must be replaced by two ghost fields, due to the non linearity of the BRST ghost field 
transformation [sc" oc f'^^'^d'c'^, see Eq. (13.81) ]. The last technical point to be clarified is the 
dependence of the STIs on the (external) momenta. One should notice that the integral over 
d'^x present in Eq. (I3.30p . together with the conservation of momentum flow of the Green's 
functions, implies that no momentum integration is left over; as a result, the STIs will be 
expressed as a sum of products of (at most two) Green's functions. 

An advantage of working with the BV formalism is the fact that the STI functional 
of Eq. (13.301) is valid in any gauge, i.e., it will not be affected when switching from one 
gauge to another. In particular, if we want to consider the BFM gauge, the only additional 
step we need to take is to implement the equations of motion for the background fields at 
the quantum level. This latter step is achieved most efficiently by extending the BRST 
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symmetry to the background gluon field, through the relations 



sA^ = n^, sn^ = 0, (3.31) 

where f2™ represents a (classical) vector field with the same quantum numbers as the gluon, 
ghost charge +1 and Fermi statistics (see also Table [T]). The dependence of the Green's 
functions on the background fields is then controlled by the modified STI functional 



where V denotes the effective action that depends on the background sources VL^ (with 
r = r'|f^=o), and >S(r')[$] is the STI functional of Eq. flCTj) . Differentiation of the STI 
functional (13.321) with respect to the background source and background or quantum fields 
will then provide the so called BQIs, which relate IPI Green's functions involving background 
fields with the ones involving quantum fields. The BQIs are particularly useful in the PT 
context, since they allow for a direct comparison between PT and BFM Green's functions. 

Finally, the background gauge invariance of the BFM effective action implies that Green's 
functions involving background fields satisfy linear Wis when contracted with the momentum 
corresponding to a background leg [see, e.g., Eq.s (12.371) - (12.401) ]. These Wis are generated 
by taking functional differentiations of the WI functional 

d'x J2 i^^^i-^)^) = 0' (3-33) 

where ■&°'{x) are the local infinitesimal parameters corresponding to the SU{3) generators 
that now play the role of the ghost field. The transformations 5^$ are thus given by 

S^c" = -gP^'c^^" 5^0" = -gr^'c^^', 

s^'iPl = igr{e),ji:i 6^iJi = -zgr^pue),,, (3.34) 

and the background transformations of the anti- fields 5^$* coincide with the gauge transfor- 
mations of the corresponding quantum gauge fields according to their specific representation. 
Notice that, in order to obtain the WI satisfied by the Green's functions involving background 
gluons A, one has to differentiate the functional (13.331) with respect to the corresponding 
parameter t?. 
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All the STIs and BQIs needed for the PT construction carried out in the rest of this 
paper, together with the method of constructing the auxiliary functions appearing in these 
identities, are reported in Appendix [D] and [El respectively. 



D. Faddeev-Popov equation(s) 

The final ingredient needed for carrying out the PT program for SDEs is the derivation 
of the so-called Faddev-Popov equation (FPE). The FPE depends crucially on the form of 
the ghost Lagrangian, which, in turn, depends on the gauge fixing function [see Eq. fl3.7p ]. 
In what follows we will first present the corresponding derivation in the gauges, and then 
in the BFM. 

To derive the FPE in the gauges, one observes that in the QCD action the only term 
proportional to the anti-ghost fields comes from the Faddeev-Popov lagrangian density, 
which can be rewritten as 



^F^G = -c'^d'^isA-) = -'^d'^j^- (3.35) 
Differentiation of the action with respect to then yields the FPE in the form of the identity 

3^7^ = 0. (3,36) 



ST 



5c^ 5A*' 



so that, taking the Fourier transform, we arrive at 

— + iq''-— = 0. 3.37 

Thus, in the i?^ case, the FPE amounts to the simple statement that the contraction of a 
leg corresponding to a gluon anti-field {A*J^) by its own momentum (g'^) converts it to an 
anti-ghost leg (c™). Functional differentiation of this identity with respect to QCD fields 
(but not background sources and fields, see below) furnishes useful identities, that will be 
used extensively in our construction. 

For obtaining FPEs for Green's functions involving BFM gluons and sources, one has to 
modify Eq. f l3.37p . in order to account for the presence of extra terms in the BFM gauge 
fixing function (and therefore in the BFM Faddeev-Popov ghost Lagrangian). Eq. fl3.35p 
gets then rewritten as 

ST' 6T'\"' 

+ '^^Tir - iP^'^^^r = 0. (3.38) 



Sc"" V SA* , 

The specific FPEs needed for the PT construction are reported in Appendix O 
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FIG. 9: The S'-matrix one-loop PT setting for constructing the gluon propagator. The external 
particles are left unspecified since they can be both quarks as well as gluons. The kernels JC^^^ 
appearing in diagram (b) are the tree-level version of the one shown in Fig.s [TT] and 1151 (depending 
on the external particles chosen) and therefore contains only the IPR terms shown there. Three 
diagrams having the kernels on the opposite side are not shown. 

E. The (one-loop) PT algorithm in the BV language 



Before entering into the intricacies of the SDEs, it is important to make contact between 
the PT algorithm and the BV formalism. This is best done at the one-loop level, since in 
this case all calculations are rather straightforward and it is relatively easy to compare the 
standard diagrammatic results with those coming from the BV formalism. This comparison 
will (i) help us identify the pieces that will be generated when applying the PT algorithm, 
and (ii) establish the rules for distributing the pieces obtained in (z) among the different 
Green's functions appearing in the calculation. 

The starting point is the embedding of the (one-loop) gluon propagator into an S'-matrix 
element (Fig. [9]), exactly as done in subsection III D[ Then, carrying out the PT decompo- 
sition F = F^ + F^ on the tree-level three-gluon vertex of diagram (b) [see Eq. f l2.12p ]. we 
find 



(b) 



(bf + {bf, 



(3.39) 



(3.40) 
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Of couse, ki and ^2 are not independent, since k2 = q — ki, thus, we have 



ram'n' ^ 
J 9afj.' 








fam'n' „ 

/ ga^l' 








fam'n' 

J iJav 









w''(^2)AX ''{h)IC\l,A^^^{ki,p2, -pi) 
'A;S't''^(^i)Ait'^(A:2)4°].^„,^(A;2,P2,-;.i), (3.41) 

and we see that the contributions of the two longitudinal momenta add up, thus removing 
the 1/2 symmetry factor (this is clearly an all-order result, since the above derivation does 
not depend on the various Green's functions and kernels being at tree- level). Therefore we 
obtain 

(bf = -gn^-g^l^ ^^k^l^%A,^^(k2,P2,-Pl)■ (3.42) 
On the other hand, using the results 

klDZ'ih) = 5"^™', ^^!^'A,Ar = 39tr''''^ rfj,^.,^^- = 0, (3.43) 



we find that the STI of Eq. flD.29p reduces to 



K^%A^^^^{h,P2,-Pi) = -gg2f--T^^^^^{p2,-Pi) + r'^'l{^^^^^^ 

+ }C%^,^^Up2,-puh)T^Z^P2). (3.44) 

At this point the calculation is over and one needs to reshuffle the pieces generated. First 

of all, notice that when the external legs are on-shell the last two terms of the above STI 
drop out, by virtue of the (all-order) equations of motion 

TfAp2Hp2) ^ = 0, (3.45) 

u{p{)T^ll{p{) = 0. (3.46) 



Thus, making use of Eq. (IE. 90 we are finally left with the result 

J ki 1 2 



-r[;i^..(-g)ri°],^(P2,-Pi). (3.47) 



Notice that, as explicitly shown in Appendix [El the auxiliary function F^^^* [see Eq. (lE.Qp 
and Fig. [2l] coincide with the function A^/? [see Eq. fl2.33p and Fig. [7] to all orders 

A,;3(g) = rn.A;(g). (3.48) 
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Thus the scalar function G{q^) introduced earher in Section Hi El corresponds also to {i times) 
the Qap part of T^^Ay 

We can now define the PT (on-shell) quark-gluon vertex, by considering the corresponding 
Green's function embedded in the diagrams 

ih? + (c) = (6) + (c) - {hf 

^ ^r«^^(p„ -p,) = .r«^^(p„ -p,) + r«^..(-g)r^°)^^fe, -p,), (3.49) 

while the PT self-energy will given by the combination 

(a) + 2{hf 11%) = n«(g) + n;:«(g). (3.50) 

The factor of 2 comes from the mirror diagram of (6) having the kernel on the left side, and 
we have defined [with the aid of Eq. (ID.lOp 

= 2<^^.(^)ri°]^.(^)- (3.51) 

We can now proceed to the comparison of the PT Green's function with that of the BFG, 
by resorting to the BQIs. Clearly, Eq. (13.491) represents the one-loop version of the BQI of 
Eq. (]E.13p . and we immediately conclude that 

r«,,-(P2,-Pi)-r«^^(P2,-Pi). (3.52) 

For the self-energy we have instead [recall that — r^m^n(A;) = n^j,(A;)] 

^"n«(^) = -r%A^^i<i) + 2^r(J)^..(g)r(°)^, (g) , (3.53) 



which represents the one-loop version of the BQI of Eq. (lE.4p . i.e., we have 



5"'nW(?) = -r^ii.(?)- (3.54) 



The procedure just described goes through almost unaltered when choosing the external 
legs of the process to be gluons. In this case 



J fci 1 2 



(3.55) 
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and one has from the STI of Eq. (1D.36|1 the result 



jj, u p a '-p-' '-a 



+ ^2.rA-A*^(P2,-pi,A;2)r5J]^„(/c2). (3.56) 



As before, the second and third terms drop out when the external gluons are taken to be 
on-shell; thus we are left with the terms 



= -r2A-(-^)ri°]A.^.(P2,-Pi) + (&')• (3.57) 



The first term is exactly the PT propagator-like piece encountered in the quark case; this is 
the essence of the process independence of the PT. Notice, however, that the second term 
was not present before. The action of this term will be discussed in detail in Section llVBt 
here it suffices to note that it is a vertex-like piece (as is evident from the structure of the 
kernel appearing in it) and, therefore, it ought to be allotted to the PT three-gluon vertex. 
Thus we can define the PT (on-shell) three-gluon vertex and propagator as before, i.e., 

(bf + {b') + (c) + (d) = (b) + (c) + (d) + (60 - {bf 

=^ ^^ALrAiiP^, -Pi) = i'^M.ArAdP^, -Pi) + rS^r (P2, -Pl), 



(a) + 2(6)P ^ n«(g) = n«(g) + Iil^^\q). (3.58) 
The comparison with the BQI of Eq. (IKTTI) shows then that 

ri'iA^Aj(P2, -Pi) = r^]^.^. (P2, -Pi), (3.59) 

and again we find 

5'''^th) = -^%X^{q)- (3-60) 

The (one-loop) procedure described above carries over practically unaltered to the cor- 
responding SDEs. This is due to the fact that: (z) the pinching momenta will be always 
determined from the tree- level decomposition of Eq. ( 12.12P : (m) their action is completely 
fixed by the structure of the STIs they trigger [Eq.s ( ID. 291) and (ID. 361) for the vertices at 
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hand]; (iii) the kernels appearing in these STIs are the same appearing in the corresponding 
BQIs; thus, it is always possible to write the result of the action of pinching momenta in 
terms of auxiliary Green's functions appearing in the BQIs. 

The only operational difference is that, in the case of the SDEs for the quark-gluon 
vertex and the three-gluon vertex, all three external legs will be off-shell. This is of course 
unavoidable, given that these (fully dressed) vertices are nested in the SDE of the off-shell 
gluon self-energy [see Figj6l diagrams (dn) and (c/i), respectively], and their legs inside the 
diagrams are irrigated by the virtual off-shell momenta. As a result, the equations of motion 
employed above [viz. Eq. (13.461) ] should not be used in this case; therefore, the corresponding 
terms, proportional to inverse self-energies, do not drop out, and form part of the resulting 
BQI. 

Thus the PT rules for the construction of SDEs may be summarized as follows: 

i. For the SDEs of vertices, with all three external legs off-shell, the pinching momenta, 
coming from the only external three-gluon vertex undergoing the decomposition (12.121) . 
generate four types of terms: one of them, corresponding to the term {b') in Eq. (13.571) . 
is a genuine vertex-like contribution that must be included in the final PT answer for 
the vertex under construction, while the remaining three-terms will form part of the 
emerging BQIs (and thus would be discarded from the PT vertex). These latter 
terms have a very characteristic structure, which facilitates their identification in the 
calculation. Specifically, one of them is always proportional to the auxiliary function 
^UA*, while the other two are proportional to the inverse propagators of the fields 
entering into the two legs that did not undergo the decomposition of (12.121) . 

ii. In the case of the new SDE for the gluon propagator the pinching momenta will only 
generate pieces proportional to Tqa*, which should be discarded from the PT answer 
for the gluon two-point function (since they are exactly those that cancel against the 
contribution coming from the corresponding vertices). 

IV. PT GREEN'S FUNCTIONS FROM SCHWINGER-DYSON EQUATIONS 

After the introduction of the useful tools and basic rules required for the application of the 
PT program to the (non-perturbative) case of SDEs, we are ready to describe in detail the 
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FIG. 10: The SDE for the quark-gluon vertex. The symmetry factors of the diagrams (first 
line) are: s{a,b) = 1. s(c) = 1/2, s{d) = 1/6 and s(e) = —1. For the key diagram (c) we show 
explicitly the kinematics chosen. In the second line we show the additional topologies present in 
the BFM version of the equation [s{f,g) = —1], generated dynamically by the PT procedure. 

actual construction, starting from the corresponding SDEs written in the Feynman gauge of 
the R^. We will first derive the new SDEs for the two vertices, F^^^ and Taaa, given that 
the calculations are easier to work out, and will then address the more complicated case of 
the SDE for the PT gluon propagator F^^. 

A. Quark-Gluon vertex 

The SDE of the quark-gluon vertex, shown in Fig. [TOl is the simplest one as far as the 
PT construction is concerned, capturing at the same time several of the essential steps that 
appear during the application of the PT algorithm to the SDEs of QCD. 

We start by carrying out the decomposition of Eq. (12.121) on the tree-level vertex appearing 
in (c), the only diagram we will touch in our construction. Let us concentrate on the F^ 
part; one has 
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FIG. 11: Skeleton expansion of the kernel appearing in the SDE for the trilinear quark-gluon 
vertex [see diagram (c) of Fig. [TU]. Black, white, and gray blobs denote IPX functions, connected 
functions, and SD kernels, respectively. 



ki 



(4.1) 



where the kernel JCaaiPiP is shown in Fig. [TTl Using the STI of the kernel /C^m^n^^ given in 
Eq. flD.291) . we obtain from (14. ip four terms, to be denoted by (si), (^2), (S3) and (54), i.e., 



(cf = (si) + (S2) + (S3) + (S4), 



(4.2) 



with 



{si) = gn'^^'gaw [ D'"''"(A;i)A:;::;(A;2)r,™^,M-(A:2,-fci-A;2)r^.^^(p2,-pi) 

J ki 

{S2) = gn'^^'gau'T^^ip,) [ D'"''"(fci)A::::;(A;2)/c^.^,™^.(P2,fci,-pi), 

Jki 

{S3) = gn'^^'go..' [ /^'"''"(A;i)A:;'4(A;2)/C^n^,^V"(P2,-Pi,fci)r^^(p2), 



(s,) = gn^'gaw / D^^{k,)^':^ak2)T,^Al^^^{h,P2,-Pi)TA^^A-Sh). (4.3) 



Next, using Eq.s flE.9p . flE.14p and flE.lSp . it is fairly straightforward to demonstrate that 



(si) = -^n-A:A-(l)'^AiM{P2,-Pi 



r^^(pi)r^Qa^. (g, —pi), 



(4.4) 



Evidently, (si) gives rise to the PT propagator-like term, while (^2) and (53) generate the 
terms that in the usual S'-matrix PT would vanish on-shell, due to the (all-order) spinor 



equations of motion T^j^{p2)u{p2) 



P2=m 



0, and u{pi)T^j,{pi) 



ii=m 



0. Of course, in 



our case we are not allowed to use the equations of motion, given that the quark legs are 
considered to be off-shell. 
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FIG. 12: The Schwinger-Dyson equation (j4.9p satisfied by the ghost propagator. 



Let us finally look at the term (54), and show how it combines with the remaining 
diagrams to generate the BFM quark-gluon vertex T^^^. To this end, using Eq. fID.llI) and 
the FPE satisfied by the IPI function F^^,^^, we write (54) = (s4a) + (54^), with 



{S4a) = -ign'^gaJ /^"'™(A:i)F,™^..^^(fc2,P2,-pi) 

J ki 



[S4b) - gn " go 



(4.5) 



'fci '^2 

The general structure of these two terms suggests that (54^) should give rise to the ghost 
quadrilinear vertex, while (54;,), when added to diagram (e), should symmetrize the trilinear 
ghost gluon coupling. It turns out that this expectation is essentially correct, but its realiza- 
tion is not immediate, mainly due to the fact that (34^) contains a tree-level instead of a full 
ghost propagator [(A^D^^ instead of D(A;2)], while (s4a) can reproduce, at most, diagram (/) 
of Fig. [ini but not (jj). The solution to this apparent mismatch is rather subtle: one must 
employ the SDK satisfied by the ghost propagator, shown in Fig [121 This SDK is common 
to both the i?g-gauge and the BFM, given that there are no background ghosts. 

To see how this works in detail, add and subtract to Eq. (14. 5 p the missing term (see 
Fig. [13]), obtaining 

(S4a) = -ign'^g., I D^'"\k,) [r,™^*.^^(fc2,P2,-Pi) 

J ki 



- T'^,A*i{k2)iD3^\k2)T^,^^,.^^{k2,P2. -pi 



-ign'^ga, [ D--\k{)ic,^,^..^^{k 

J ki 



■2,P2,-Pl; 



lS46j 



-gn gl 



dn' k2v' 

h.2 " caA 



X L'™'™(A;i)r,™gd^^(A;2,p2,-pi 



(4.6) 



(4.7) 



where the auxiliary function F^^* has been defined in Eq. f lE.7p . and is given by Tca* minus 
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FIG. 13: Diagrammatic decomposition of the kernel /C^m^*7^^. The second term is the one added 
(and subtracted) to the original sum (s^a) + (sib) of Eq. (j4.5p . After replacing the gluon anti-field 
by the corresponding composite operator (second line) , taking into account the extra structure 
provided by the (c)^ term, this kernel furnishes the BFM terms (/) + (g) [see Eq. (14. 8p .] 



its tree-level part. Using Eq. flE.Sp . we can then rewrite {s^a) as 



'ki Jki 



X 



A''' /C'^' 
if) + (?)• 



(/C3, h,P2, -Pi) - ir^3^.'ge'(A;3, h)D'^'^'{k2)T^,^^,'^ T{k2,p2, -pi 



(4.8) 



We next turn to (^4^) and consider the ghost SD equation of Fig. [121 One has 

xdn' Sdg 



(4.9) 



where, as before, T'^^.^, is given by T^g^gi minus its tree-level part. Multiplying the above 



equation by fcg, using the FPE (1C.1|) and factoring out a k2u' we get the relation 



k2u'D''^'{k2) = 6" 



K.A^Ak2)D^-\k2) 



6' 



'^2 

'''''j^-K.A*Ak2)D^''ik2). 

run u' 



(4.10) 



Therefore, we obtain 



{sAb) = -9n'''' f k2aD"^'"^ik^)D-'-ik2)T,^,„^^{k2,P2,-pl). (4.11) 
Jki 
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Adding this last contribution to diagram (e) we finally arrive at 

(e) + (s4fe) = ^^r"'"' [ {h-h)aD'^''^{h)D^'^{k2)IC„4h,p2,-pi 



= (e). (4.12) 

Next, observe that the graphs (a), (6), and (d) of Fig. [10] can be converted to hatted 
ones automatically (see the corresponding tree-level Feynman in Appendix [Fj, and that 
(c)F = (c) since in the BFG = vf^^. Thus, 

- r^^(pi)F^^.^.(g, -p,) + [(a) + (6) + (c) + (d) + (e) + (/) + (?)]^. 

(4.13) 

The sum of diagrams in the brackets is nothing but the kernel expansion of the SDK gov- 
erning the vertex F^^^, i.e., 

^ri.^^(p2, -Pi) = [(a) + (b) + (c) + id) + (e) + (7) + mi (4.14) 

After this identification, it is clear that Eq. (I4.13P coincides with the full BQI of Eq. ( lE.lSp . 
namely 

In summary, the application of the FT to the conventional SDK for the quark-gluon 
vertex {i) has converted the initial kernel expansion [graphs (a) to (e) in Fig. [10] into the 
graphs corresponding to the kernel expansion of the vertex F^^^; (ii) all other pinching 
terms extracted from the original diagram (c) are precisely the combinations of auxiliary 
Green's functions appearing in the BQI that relates the initial vertex F^^^ with the final 
vertex F^^^. 

Notice at this point that the skeleton expansion of the multi-particle kernels appearing 
in the SDE for F^^^ is still written in terms of the conventional fully dressed vertices and 
propagators (involving only quantum fields). Thus, Eq. (14.141) is not manifestly dynamical, 
i.e., it does not involve the same unknown quantities on the right and left hand side; this 
situation is exactly analogous to the gluon propagator case discussed in subsection III El 
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Specifically, in order to convert (14.141) into a genuine SDE, one has two possibilities, both 
involving the use of the above BQI: (z) substitute the Ihs of Eq. (14.151) into the rhs of 
Eq. (14. 141) and solve for the conventional F^^^ vertex, or (ii) invert Eq. (14.151) and use it to 
convert every F^^^ vertex appearing in the rhs of Eq. (I4.14p into a F^^^ vertex. It would 
seem that the latter option is operationally more cumbersome, especially taking into account 
that a similar procedure has to be followed for all the Green's functions that appear in the 
coupled system of SDEs that one considers. 



B. Three-gluon vertex 

The construction of the FT three-gluon vertex proceeds in a very similar way, with some 
additional subtleties that we will spell out in detail in what follows. We emphasize that 
the purpose of this exercise is to generate dynamically the vertex F^^^ and not the fully 
Bose-symmetric vertex F^^^ studied in 15|, |29|]. The reason is that it is the former vertex 



that appears in the SDEs for the gluon propagator [see, e.g., diagram (di) in Fig. [6], making 
it the relevant object to consider at this level. 

We start by considering the conventional SDE for the three gluon vertex (Fig. [HI), 
and carry out the standard F^ + F^ decomposition to the tree-level vertex of diagram (c) , 
which is the only one we will modify in our construction. 

We then find 

{cf = -^^r™'"' / (wW -^a/.'^2i.0AmCi(^l)^n'n(^2)/CA™AlJA-A|(A;2,P2,-Pl) 



2^ 



1 



= gn'^'go.u' I ^^irMK^A-A-^ArAsSk2.P2.-vi). (4.16) 

J k\ 1 

and the kernel Kaaaa is shown in Fig. [T51 

The next step is to apply the STI of Eq. flD.361) . and scrutinize the various terms, denoted 
again by (si), (^2), (S3), and (54). 

For the first three terms, we get the following results 

(^1) = gn'^-'g^w I i)'"''"(A;i)A::::;(fc2)F,™^,^*.(fc2,-fci-fc2)F^.^.^.(P2,-Pi) 

J ki 

= -rngA7(-9)rAdA-Aj(p2, -pi), 
(■52) = gT'^'^^'gau' / ^™'™(/i;i)Aj;'^(/i;2)^c™A;)AjA7(^2, -pi,P2)FAdA-(p2) 

J ki 

= -rn.A.A7(-Pl,P2)F^d^.(p2), 
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FIG. 14: The SDE of the three-gluon vertex. The symmetry factors of the (first and second 
hne of the figure) diagrams are s{a,b) = 1. s{c) = 1/2, s{d) = 1/6, s{e) = —1, s{f,g) = 1/2. 
In the third line we show the additional topologies present in the BFM version of the equation 
[s{h,i,l,'m) = —1], generated during the PT procedure. 

(ss) = gn'^-'g^w f /^™''"(A;i)A-:(fc2)/C,™AijA;;A7(^2,P2,-Pi)r^.Aj(Pi) 

= -rngA-A7(p2,-Pi)r^dAj(pi)- (4.17) 

As in the case of the quark-gluon vertex, (si) represents the propagator-like contribution 
that in the S'-matrix PT would be allotted to the new two point function F^^. Notice that 
this term is equal (modulo the external vertex) to the one found in the quark-gluon vertex 
case; this is the manifestation of the well-known property of the process-independence of 
the PT algorithm (as already noticed in our previous one- loop analysis): the propagator-like 
contributions do not depend on the details of the external (embedding) particles. As for (52) 
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FIG. 15: Skeleton expansion of the kernel appearing in the SDE for the three- gluon vertex [diagram 
(c) of Fig. [H] 

and (ss), they correspond again to terms that would vanish on-shell, but now are retained 
in the final answer due o the off-shell condition of the external legs. 
Finally, one has to consider the term (54), given by 

(54) = gn'^^'g^w [ /^™'"(fci)A:;::;(A:2)/c,™^*.^.^.(P2,-Pi,A:2)r^.A.>(fc2), (4.18) 

J ki 

which again can be written as the sum of the following two terms 

{s,a) = -ign'^g., [ /^™''"(fci)/C,.^..^.^.(fc2,P2,-Pi), (4.19) 

J ki 

{s,b) = -gn'^^'gaw I 5'-'^D^'^{h)}C,,.,.^rAi{k2,P2,-Pi). (4.20) 

hi 2 

The kernel JCccAA is defined by replacing in Eq. flD.39p every anti-field leg A* by the corre- 
sponding anti-ghost field c. As before, (s4b) has a tree-level ghost propagator, while (s^a) 
misses a diagram that we need to add and subtract to solve the two problems simultaneously. 
Even so, we are still missing the diagrams (/) and (m) of Fig. [TH they will be generated by 
the tree-level contribution appearing in the SDE of the auxiliary function T^aa*- In order 
to isolate this contribution as early as possible, let us write 

^c^A:''ArAs {k2,P2,-Pl) = (A;2,P2,-Pi) 

+ iT,^^,,A-pij)^D''m'^XrA-'^p^^^^) 

= K^A'''ArAsih,P2,-Pl)-^gf''''g;r^^^s,,{-p,,i)D^'\i) 
- ^gf''g2D^'\nrc-Ar-cAP2,-n, (4.2i) 

where the prime denotes that the TcAA* that appears in the IPR terms starts at one-loop. 
We then find (see also Fig. [16]) 

(S4a) = -^gn'''ga, [ /^™''"(A:i)/C^.^..^.^.(fc2,P2,-pi) 

Jkl d p ^ 
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FIG. 16: The IPR terms appearing in A^c'"A*"'A^A'' contain a tree-level contribution generating the 
missing BFM topologies. Here we show the case for (/); a symmetric term generates (m). The first 
term on the rhs is part of the skeleton expansion of diagram (h) of Fig. [TH Notice that the Ihs is 
simply a pictorial representation of the rhs, taking advantage of the notation introduced in Fig|13t 
the anti-fields are static sources and do not propagate. 



-Pl, 



Jki 

Jki ^ 

= (4a) + (0 + (^)- (4.22) 
For generating the remaining terms one proceeds as in the quark case, writing (see Fig. [T7|) 



Jki I d P ^ 

- r^<7A7(^2)«^^^'(/i;2)/C^".g<,'A-Aj(^2,P2, -pi) 

= -ign'^g., I D"^'"'ih)ic'^\.,^^^^ik2,P2,-Pi) 

J k P (7 



1^1 ^,1 



M = -gn'^gl 



/ k2u' 



X 



X L)"'™(A;i)/Cc™gd^.4j(A;2,P2,-pi). 



(4.23) 



(4.24) 



Then, using Eq. flE.SP (which can be safely done now, since tree-level contribution has been 
already taken into account) and Eq. (14.101) . one finds 



so that 



{S4b) = -gn^ / k2c.D"'"'{k,)D^^ik2)IC,^,^Ar^A^^ik2,P2,-pi), 



{S4b) + (e) = (e). 
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(4.25) 
(4.26) 

(4.27) 




FIG. 17: Diagrammatic decomposition of the SD kernel K.^^.*^ defined in Eq. (14.230 . The first 
term represents the kernel fC' m A*-r Ar as ^ therefore the TcAA* appearing in the corresponding IPR 
terms start at one-loop. The second term is the one added (and subtracted) to the original sum 
(s4jj) + (54;,). After replacing the gluon anti-field J^^ with the corresponding composite operator 
(second line), this kernel generates the BFM terms (K) + {%). 

Using the tree-level Feynman rules (see Appendix [F]), it is straightforward to establish 
that the graphs (6), (rf), (/), and {g) can be converted to hatted ones automatically, and that 
[cf = (c). Thus, collecting all the pieces we have, and using the standard FT decomposition 
f l2.12p on the tree- level contribution (a), we get 

iTA-ArA^^{p2,-pi) = -rnaA7(-9)rAdA-Aj(P2,-pi) " AjA^ (-Pl, P2)rAdA- (^2) 

- ^niArA'/iP2, -Pi)rAdAj(Pi) + [(a) + (b) + (c) + (cT) + (e) 

+ (/) + (?) + (h) + (T) + (f) + (m)]-; - tgr'T^{p2, -pi). (4.28) 

As in the previous case, the sum of diagrams in the brackets is nothing but the kernel 
expansion of the SDK governing the vertex T^^^, i.e., 

i^MA^A^ (P2, -Pi) = [{a) + (b) + (c) + (d) + (e) 

+ (/) + (?) + (h) + (T) + (0 + (m)]-;. (4.29) 

This in turn implies that Eq. fl4.28p represents the BQI of Eq. flE.171) up to the last (tree- 
level) term in the rhs. Of course this tree-level discrepancy is to be expected since the FT 
algorithm cannot possibly work at tree-level if the external legs are amputated, as is the 
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case in the SDEs we are considering. To be sure, if we start from the tree-level T^^j^^ only, 
i.e., without hooking (two of) the external legs to (conserved) external currents, we can still 
carry out the decomposition of Eq. (12.121) . but the term will have nothing to act upon. 

Notice finally that the discussion following Eq. (14.151) the SDE for the gluon-quark vertex 
applies with minimal modification to the three-gluon vertex case discussed here. 

C. The gluon propagator 

In this section we turn to the SDE of the gluon self-energy. From the technical point 
of view the construction is somewhat more involved compared to that presented for the 
vertices, simply because the FT decomposition of Eq. (I2.12p must be carried out on both 
sides of the self-energy diagram. Fut in a different way, now we must convert to background 
gluons not one but two external gluons. To the best of our knowledge, the most efficient 
procedure to follow consists of the three basic steps described below |l3| . 

1. First step 

The starting point is diagram (ai) of FigO Following the FT procedure, we decompose 
the tree-level three-gluon vertex according to (12.121) and concentrate on the pinch part. We 
then get 

= ign'^'g..' ^^A::ak2)KT^^^„^.^{k,,-q). (4.30) 

At this point the application of the STI of Eq. (ID. 91) together with Eq. ( ID. Ill) and the 
FFE (IC.Sp . results in the following terms 

J ki 

+ gn'^g., [ z^™'™(fci)F,.^.^..(-g,fc2) 

J ki 

Jk, M p 

= [si] + {S2) + (S3). (4.31) 
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FIG. 18: Schwinger-Dyson equation satisfied by the gluon self-energy — r^^. The symmetry factors 
of the diagrams are 621 ^e) = 1/2, s{b^) = 1/6, and all the remaining diagrams have s = —1. 



Clearly, using the SDE of the auxiliary function Fj^^. , shown in Eq. flE.91) . one has immedi- 
ately that 

(si) = -iFj^g^*7(g)F^d^6(g). (4.32) 

This would be half of the pinching contribution coming from the vertex in the S'-matrix FT. 

As far as the (S2) and (53) terms are concerned, let us start by adding and subtracting 
to them the expression needed to convert the tree-level ghost propagator of (53) into a full 
one; making use of the ghost SDE fl4.10p we then get 

{S2) = -gn'^Qa, [ zD"^'"^{h) [^F,„^.^..(-g,fc2) + r,,^..(A;2)/^^'^(fc2)r,™A^g.(-g,fc2) 

Jki I 13 d a fi 



^3 



J ki 



The second term symmetrizes the trilinear ghost-gluon coupling, and one has 



(S3) + (03) = ih), 



(4.33) 



(4.34) 



where (63) is shown in Fig. [181 The term (52) will finally generate all the remaining terms. 
To see how this happens, we denote by (s2a) and {s2b) the two terms appearing in the square 
brackets of (s2a), and concentrate on the first one. Making use of the SDE (lE.Sp satisfied 
by the auxiliary function Tcaa* and the decomposition ( IE. Ill) of the kernel appearing in the 
latter, we get 



{S2a) = g^n'T^'g^p [ D"^'"^{h 

Jki 
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FIG. 19: Schwinger-Dyson equation satisfied by the gluon self-energy — r^^. The symmetry factors 
of the diagrams are s(ci,C2,C6) = 1/2, 5(03,04,07) = —1, 5(05) = 1/6. 

= (64) + (67) + %) + (fell). (4.35) 
Using instead the SDE satisfied by Fca*, shown in Eq. (IE. 70 . we obtain 



Xrcm^tcsl"?) ^2) 



(&9 



(4.36) 



Finally, since diagrams (02), (04) (0.5) and (ag) carry over to the corresponding BFM ones 
(^2), (^5), (^e) and (611) and {a\f = (bi), we have the final identity 



(S2) + (S3) + 



1=2 



11 



EC" 



1=1 



and therefore 



^A?.A''S^) = -^rs^a^*7(g)F^d^b (g) - F;j„^6 (g), 



7^/3 



(4.37) 



(4.38) 



which is the BQI of Eq. (jM]). 



2. Second step 



The second step in the propagator construction is to employ the obvious relation 



(4.39) 
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that is to interchange the background and quantum legs (the SDE for the self-energy — is 
shown in Fig. [T9l) . This apparently trivial operation introduces a considerable simplification. 
First of all, it allows for the identification of the pinching momenta from the usual PT 
decomposition of the (tree-level) F appearing in diagram (ci) of Fig. [12] [something not 
directly possible from diagram (&i)]; thus, from the operational point of view, we remain on 
familiar ground. In addition, it avoids the need to employ the (formidably complicated) BQI 
for the four-gluon vertex; indeed, the equality between diagrams (cs), (ce), (cy) of Fig. [19] 
and ((is), (rfg), (c^n) of Fig. [6], respectively, is now immediate [as it was before, between the 
diagrams (04), (05), (ag) and (65), (fee), (^ii), respectively]. 

3. Third step 

We now turn to diagram (ci) and concentrate on its pinching part, given by 

{cir = ign'''gau' [ ^,A';:ak2)KT^^^„^Jh,-q). (4.40) 

Jki '^l M /J 

Notice the appearance of the full BFM vertex F^^^ instead of the standard Taaa (in the 
R^). The STI satisfied by the former vertex has been derived in Eq. flD.25p . Now, the first 
three terms, (si), (^2) and (ss), appearing in this STI, will give rise to PT contributions 
exactly equal to those encountered in first step described above, the only difference being 
that the field appearing there is now a background field A^. Thus, following exactly the 
reasoning described before, we find [see again Fig. [6] for the diagrams corresponding to each 

{d^)] 

(si) ^ -iFn.A-(g)r^ei6(g), (4.41) 

is2) + (S3) + (C3) = (d^) + (d,) + (dj) + (ds) + (dg) + (rfio). (4.42) 
For the term (54) we have instead 

is,) - g'r'^'^r'-'ga.'g^, [ A^^Liki). (4.43) 

Jki 

Clearly this has a tadpole-like structure; in particular, it is immediate to prove that when 
added to (02) it will convert it into (^2) 

(S4) + (C2) = (^2). (4.44) 
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FIG. 20: Summary of the PT procedure employed in the text in order to construct the new PT 
SDE of the gluon propagator. 



Thus, since as always (ci)^ = (di) we get 

{S2) + (S3) + (54) + 

and therefore 



i=2 



11 



(4.45) 



i=l 



- rAsA^(?) = -^rnsAr(?)rAeA^(g) - r^.^^ (g), (4.46) 

which is the BQI of Eq. (lE.Sp . This concludes our proof. 

In Fig. [20]we summarize the steps that allowed the successful construction of the SDE for 
the PT propagator; putting together the three steps above, we have been able to generate 
the complete BQI of Eq. (IE. 41) . namely 

^r^gib (g) = iT^aAo^iq) + 2Ff^.^*7(g)F^d^6 (g) - irj^„^.7(g)F^d^.(g)Ff^6^«(g). (4.47) 

According to the PT rules put forward in Section IIII El on the one hand the PT gluon 
two-point function zTaa would coincide with the rhs of Eq. f l4.47p after dropping the terms 
proportional to the auxiliary function Tqa*, since these would cancel anyway after adding 
the contribution coming from the corresponding vertices. On the other hand, recalling that 
^n^A* coincide with A^^, and observing that the only relevant part in the identity above 
of such functions is the one proportional to the metric tensor (due to the transversality of 
the gluon two-point function Taa), it is immediate to show, using Eq. (12.61) and the relation 



— r^^^^ = Ila(3, that Eq. f l4.47p can be cast in the form of the SDE shown in Eq. fl2.36p . 



D. How to truncate the new SDEs 

After constructing the new SDE series, let us focus on its truncation properties, both 
from the theoretical as well as the practical point of view. 
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As has been stated repeatedly, the main theoretical advantage of the new SD series is that 
it allows for a systematic gauge-invariant truncation, in the sense described in subsection III E[ 
There we focused on how to truncate the SDE for the gluon self-energy, shown in FigEl by 
exploiting the fact that the various fully dressed graphs organize themselves into gauge- 
invariant subsets [those appearing in Eq. fl2.41l) ]. The practical importance of this property 
is the following: one can reduce the number of coupled SDE that one must include in order 
to maintain the gauge (or BRST) symmetry of the theory intact, as reflected, for example in 
the validity of Eg. (11. II) . Thus, in the case of pure Yang-Mills, within this new formulation, 
the minimum number of equations that one must consider is only two: The SDE for the gluon 
self-energy, given by the first gauge-invariant subset only {i.e., [{di) + {d2)]^f^ in FigE]) and 
the SDE for the full three-gluon vertex, shown in Fig{T3| (which is instrumental in assuring 
the gauge invariance of the subset chosen). This is to be contrasted to what happens within 
the conventional formulation: there the SDEs for all vertices must be considered, or else 
Eq. ( 11. ip is violated (which is what usually happens). 

Notice an important point, however: the present analysis does not furnish a simple dia- 
grammatic truncation, analogous to that of the gluon self-energy, for the SDE of the three 
gluon vertex ^ ^ {ki,k2), shown in Fig{TH Thus, if one were to truncate the SDE for 
the three-gluon vertex by keeping any subset of the graphs appearing in FigJT^ one would 
violate the validity of the all-order Wl of Eq. (I2.37P : this, in turn, would lead immediately 
to the violation of Eq. (11.11) . thus making the entire truncation scheme collapse. 

The strategy one should adopt is instead the following. Given that the proposed trun- 
cation scheme hinges crucially on the validity of Eq. (12.371) . one should start out with an 
approximation that manifestly preserves it. The way to enforce this, familiar to the SDE 
practitioners already from the time of QED, is to resort to the "gauge-technique" j47|. 
namely "solve" the WI of Eq. (12.371) . Specifically, one must express the three-gluon vertex 
as a functional of the corresponding self-energies, in such a way that (by construction) its 
WI is automatically satisfied. For example, an Ansatz with this property would be 

^a^a^aA^, h) = rf^^^Jh, h) - ^^^^ - U^Uh)] ; (4.48) 

contracting the rhs with qa = (^1 + ^2)0 yields automatically the WI of Eq. (I2.37p . Thus, the 
minimum amount of ingredients for initiating a self- consistent non-perturbative treatment 
is the SD for the gluon self-energy, consisting of [{di) + {d2)]^p, supplemented by an Ansatz 
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for the three-gluon vertex like the one given in (14.481) . Note that the "gauge-technique" 
leaves the transverse {i.e., automatically conserved) part of the vertex undetermined. This 
is where the SDE for the vertex enters; it is used precisely to determine the transverse 



parts. Specifically, following standard techniques [29|, |3l| , one must expand the vertex into 
a suitable tensorial basis, consisting of fourteen independent tensors, and then isolate the 
transverse subset. This procedure will lead to a large number of coupled integral equa- 
tions, one for each of the form-factors multiplying the corresponding tensorial structures, 
which may or may not be tractable. However, at this point, one may simplify the resulting 
equations {e.g., linearize, etc) without jeopardizing the transversality of Hfj^u, which only 
depends on the "longitudinal" part of the vertex, i.e., the one determined by fl4.48l) . Thus, 
the transverse parts will be approximately determined, but gauge invariance, as captured 
by q^H^u = 0, will remain exact. 

Note by the way that the methodology described above constitutes, even to date, the 
standard procedure even in the context of QED, where the structure of the SDE is much 
simpler, given that the SDE for the photon contains one single graph [diagram (ae) in FigH], 
and the photon-electron vertex satisfies automatically a naive all-order WI. Thus, while the 
PT approach described here replicates QED-like properties at the level of the SDEs of QCD, 
in our opinion a striking fact in itself, does not make QCD easier to solve than QED. 

The reader should appreciate one additional point: any attempt to apply the approach 
described above in the context of the conventional SDE is bound to lead to the violation 
of the transversality of H^u, because {i) the vertices satisfy complicated STI's instead of 
the Wis of Eq. f l2.37l) - fl2.40l) . a fact that makes the application of the "gauge-technique" 
impractical, and {ii) even if one came up with the analogue of Eq. (14.481) for all vertices, 
one should still keep all self-energy diagrams in Fig{T]to guarantee that q'^Hf^u = 0. From 
this point of view, the improvement of the present approach over the standard formulation 
becomes evident. 

Finally, one should be aware of the fact that there is no a-priori guarantee that the gauge- 
invariant subset kept {i.e., [{di) + {d2)]a/3) capture necessarily most of the dynamics, or, in 
other words, that they represent the numerically dominant contributions (however, for a 
variety of cases it seems to be true 4^). But, the point is that one can systematically improve 
the picture by including more terms, without worrying that the initial approximation is 
plagued with artifacts, originating from the violation of the gauge invariance or of the BRST 
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symmetry. 



V. CONCLUSIONS 

In this article we have presented a detailed derivation of a new SD series for non-Abelian 
gauge theories, based on the PT and its correspondence with the BFM. The procedure 
we followed for constructing the PT SDE is identical to that followed in the perturbative 
construction, without any additional new assumption. 

Our starting point is the conventional SDEs for the vertices and the gluon self-energy 
written in the Feynman gauge. The first step in the derivation is to simply carry out the PT 
decomposition of the elementary (tree-level) three-gluon vertex, Tafj,u given in Eq. fl2.12p . to 
the external three-gluon vertices appearing in the corresponding diagrams of the standard 
SD series. The part of Ta^^ denoted by F^^^, contains longitudinal momenta which get con- 
tracted with the kernels or the fully-dressed Green's functions appearing inside the diagram 
containing the original Tafiu, triggering the corresponding STIs. These STIs, in turn, contain 
pieces that, according to the well-established PT criteria, either form part of the answer, in 
this case the diagrammatic expansion of the SDE for the corresponding Green's function, 
or they are discarded from it. In Section [IV] we have worked out in detail three cases: (i) 
the quark-gluon vertex, {ii) the three-gluon vertex, and [iii) the technically more involved 
case of the gluon self-energy. It turns out that the diagrams comprising the PT answer are 
identical to those corresponding to the BEG. Thus, the resulting new SDEs, generated after 
the characteristic PT rearrangements have taken place, and the PT criterion for identifying 
the answer has been employed, correspond to the BFM SDEs, written in the BFG. This 
is an important result, because it proves the PT-BFG correspondence at the level of the 
SDE of the theory; obviously, all results on this point presented in the literature so far are 
included in the result presented here, given that any order in perturbation theory is already 
contained in the SDEs we consider. 

An additional important result, obtained from the same procedure, is the diagrammatic 
derivation of the BQIs, which, to date, have only been formally derived in the context of the 
BV formalism. The way that the terms comprising the BQIs appear in the present analysis 
is automatic: they are simply the leftovers of the PT construction, i.e., the pieces that have 
been discarded from the PT answer. 
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As explained in detail in 13|, |37|] , and mentioned also in Section [TTl the new SDE series 



for the gluon self-energy contains fully dressed vertices that satisfy simple, QED-like [i.e., 
tree-level-like) Wis, instead of STIs. This fact allows for the truncation of the SDE series 
while maintaining the transversality of the answer at any step. 

Returning to the construction of the gluon SDE, a crucial ingredient for the proof has been 
the interchange of background and quantum legs done in subsection IIV C 21 i.e., T^^ T^-^ 
(see also Fig. [201) . As we have mentioned there this allowed us to [i) identify the pinching 
momenta from the usual PT decomposition, and (ii) avoid the use of some otherwise indis- 
pensable multi-leg BQIs. But more importantly it unveils a recursive pattern that can be 
used to generalize the construction to n-point Green's functions with n arbitrary. Work in 
this direction is already in progress. 

As we have emphasized in the Introduction, the BFG is singled out dynamically when 
carrying out the PT rearrangement of a physical quantity, such as an S'-matrix element 
or a Wilson loop. In particular, after the full cancellation of all (effectively propagator- 
like) gauge-dependent pieces has taken place, and after the vertices have been forced to 
obey Abelian Wis, the resulting self-energy contribution (to be identified with the PT self- 
energy) coincides with the BFM gluon self-energy, calculated in the BFG. In that sense 
the BFG is very special, because it captures the net gauge- independent and universal {i.e., 
process-independent) contribution contained in any physical quantity. In practice, however, 
one would like to be able to truncate gauge-invariantly {i.e., maintaining transversality) sets 
of SDEs written in different gauges. This becomes particularly relevant, for example, when 
one attempts to compare SDE predictions with lattice simulations, carried out usually in 
the Landau gauge. One of the most powerful features of this formalism, not explored in 
this article, is that it can be generalized to any other gauge choice. In particular, Eq. fl2.36p 
maintains the same form, regardless of the gauge chosen. The way to accomplish this is to 
use the "generalized" PT, developed in 4l|]. The generalized PT modifies the starting point 
of the PT algorithm, namely Eq. fl2.12p . distributing differently the longitudinal momenta 
between F^^^ and F^^^. Specifically, the non-pinching part, i.e., the analogue of F^^^,, 
must satisfy, instead of (12.131) . a WI whose rhs is the difference of two inverse tree- level 
propagators in the gauge one wishes to consider. The way this works is the following. One 
starts out with the conventional SDE in the chosen gauge, carrying out the generalized 
PT vertex decomposition. Then, the action of the corresponding F^^^ projects one to the 
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corresponding BFM gauge; this includes the covariant gauges, such as the Landau gauge, 
or even non-covariant gauges, such as axial or light-cone gauges (for the way how to use 



non-covariant gauges within the BFM framework, see j41|). This new SD series contains 
full vertices that, even though they are in a different gauge, satisfy the QED-like Wis given 
in Eq.s (12.371) - (12.411) . Therefore, the truncation properties of this SDE are the same as 
those discussed in Section |TT] for the case of the Feynman gauge. The analogy is completed 
by realizing that the BQIs in the corresponding gauge allow one to switch back and forth 
from the conventional to the BFM Green's function. Thus, one may obtain, for example, 
transverse approximations for the gluon propagator in the conventional Landau gauge by 
studying the SDE written in the BFM Landau gauge, computing the [1 + G'(g^)]^ in the 
same gauge, i.e., employ Eq. (I2.36P using for the diagrams on its rhs the BFM Feynman 
rules in the Landau gauge (see Appendix). This Landau gauge SDE has already been used 



m [42] , in order to derive results for the gluon and ghost propagators that are in qualitative 
agreement with recent lattice data; as explained there, particular care is needed when taking 
the limit C,q 0. 

It would be important, both from the theoretical point of view as well as for the prac- 
tical applications, to study in detail the renormalization properties of the new SDEs (for a 
general discussion see the Appendix [B]) . For addressing this problem the intrinsically non- 
perturbative renormalization method known as "displacement operator formalism" [43] may 
prove to be particularly suitable. 

It is also very appealing to believe that the SDE derived here may be actually obtained 
from a variational principle, i.e., as a result of the extremization of an appropriate effective 
action, as happens in the case of the CJT formalism [3]. Calculations in this directions are 
already in progress. 

Finally, it would be interesting to establish connections between the dynamics obtained 
from the SDEs derived here and results based on the non-perturbative BFM formalism 
developed in [4^. 
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APPENDIX A: WHICH WAY TO PINCH AND WHY 



Historically, the basic conceptual difficulty associated with the generalization of the PT 
beyond one loop has been to determine the origin of the pinching momenta. Let us assume 
that, without loss of generality, one chooses from the beginning the conventional Feynman 
gauge. Then, the only sources of possible pinching momenta are the three-gluon vertices. 
The question is whether all such vertices must be somehow forced to pinch, or, in other words, 
whether the standard PT decomposition of Eq. (12.121) should be carried out to all available 
three-gluon vertices. The problem with such an operation, however, is the following: for the 
case of a three-gluon vertex nested inside a Feynman diagram, how does one choose which is 
the "special" momentum? Or, in other words, which way is one supposed to brake the Bose- 
symmetry of the vertex? Turns out that the solution to these questions is very simple 30|: 



one should apply Eq. (12.121) only to the vertices that have the physical momentum incoming 
(or outgoing) in one of their legs (not mixed with virtual momenta); the special leg is 
precisely the one carrying q, i.e., the physical momentum transfer appearing in the problem. 
We will call such a vertex "external". All other vertices are not to be touched, i.e., they 
should not be decomposed in any way; such vertices have virtual momenta entering into 
every one of their three legs, and are called "internal" (see Fig j2T|) . 

The reason why all other three-gluon vertices inside the loops should remain unchanged 



(no s 
PT 211 



ittin g) can be best understood by resorting to the absorptive construction of the 



22 



231]. The basic philosophy behind the absorptive construction is to emulate as 



much as possible the text-book reconstruction of the real part of the vacuum polarization of 
QED (containing say a muon-loop) from the tree-level cross-section for e"*"e^ —>■ , i.e., 
the optical theorem, and a (once-subtracted) dispersion relation. As in QED, in the case 
of the PT the basic observation also happens already at one loop: the PT subamplitudes 
(self-energies, vertices, boxes) satisfy the optical theorem individually, in a way similar to 
what happens with scalar theories and QED. 

Specifically, let us write S" = 1 -|- iT, and consider the forward scattering process 
(1{Pi)q{P2) <liPi)<l{P2), with s = q"^ = {pi + ^2)^- Restricting ourselves to only gluonic 
intermediate states, the PT amplitudes, at lowest order satisfy 

53m(gg|TW|gg), = ^ x i / [{qq\T^%g){9g\T^%qr]i, (A.l) 
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FIG. 21: Left panel: Some examples of external ( indicated with an arrow) and internal vertices 
(all the remaining). Diagram (a) has only internal three-gluon vertices, while diagram (b) has two 
internal vertices and one external, indicated by the arrow. Right panel: The two- and three-particle 
Cutkosky cuts (cutting through gluons only). 

with Jpg^ denoting the two-body phase space for massless gluons. In the equation above the 
superscript [n] denotes the order of the corresponding amplitude in the coupling constant 
g (when counting powers of g remember the couplings coming from the vertices with the 
external particles); the subscript i = 1,2,3 denotes, respectively, the propagator-, vertex-, 
and box-like parts of either side (to recover the full optical theorem, one simply sums both 
sides over i). Finally, the extra factor of | is statistical, since the final state gluons are 
considered as identical particles in the total rate. 

The meaning of propagator-, vertex-, and box-like is clear as far as the Ihs of Eq. (lA.ip 
is concerned: one must determine the imaginary (absorptive) parts of the three one-loop 
FT subamplitudes obtained after casting the amplitude q{pi)q{p2) (i{Pi)q{P2) into the 
FT form, following the standard FT rules. To get these absorptive parts one may carry 
the corresponding Cutkosky cuts to the various integrals, (including "unphysical" contribu- 
tions coming from ghost loops) or, equivalently, study where the various logarithms develop 
imaginary parts. 

Let us see now what propagator-, vertex-, and box-like means on the rhs, consisting of 
the squared amplitude for the tree-level process q{pi)q{p2) — ^ g{ki)g{k2) (with ki and ^2 
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integrated over all available phase space). A PT-rearranged squared amplitude means the 
following. Consider a normal squared amplitude, i.e., the product of two regular ampli- 
tudes [remember that "product" means that they are also connected (multiplied) by the 
corresponding polarization tensors]. Then each amplitude must be first cast into its PT 
form, by again simply following the standard PT rules. However, this is not the end of the 
story as far as the PT-rearrangement of the square is concerned. One must go through the 
additional exercise of letting the longitudinal momenta coming from the polarization vectors 
trigger a particular cancellation between the s-channel and the t-channel graphs (known in 
the literature as the "s-t cancellation"). That will finally identify the genuine propagator-, 
vertex-, and box-like pieces of the entire product. 

Now the important step is the following: Suppose that one starts out with the rhs of 
Eq. flA.l|) . i.e., one works at the level of the physical squared amplitude. The PT rear- 
rangement of the rhs may furnish the PT rearranged amplitudes on the Ihs, through an 
(appropriately subtracted) dispersion relation. Thus, the absorptive PT construction means 
to (z) PT-rearrange the rhs, {ii) impose the optical theorem (individually for each i), and 
(Hi) use analyticity to get the real parts of the PT amplitudes. 

Let us now see how the PT absorptive construction gets generalized to higher orders, 
and, in particular, how it can furnish a unique way for defining the PT construction without 
any a-priori reference to the BFM and its special vertices. 

At the next order in Eq. (1A.1I) becomes 

Qm{qq\T^%q)e = ^ (^) 2^e[{gg\T^\q)* {gg\T^\q)], 

+ \ (^) [{999\T^%q)*{999\T^%q)]i. (A.2) 

where now Jpg^ denotes the three-body phase-space for massless gluons. 

According to Eq. (]A.2|) then the imaginary parts of the two-loop PT Green's functions 
(under construction) are related by the optical theorem to precisely identifiable and very 
special parts of the squared amplitudes for the processes qq — gg and qq ggg. In 
particular, the two-particle Cutkosky cuts of the two-loop PT self-energy are related to the 
propagator-like part of the PT-rearranged one-loop squared amplitude for qq —>■ gg, while, 
at the same time, the three-particle Cutkosky cuts of the same quantity are related to the 
propagator-like part of the PT-rearranged tree-level squared amplitude for qq —>■ ggg. The 
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FIG. 22: The product of the PT-rearranged amphtudes of the process qq — > qq at one-loop (up) and 



tree- level (down). The longitudinal momenta from the 
cancellations between s-channel and t-channel graphs 
rhs of Eq. (|A2]). 



polarization tensors will produce additional 



3, 



furnishing finally the first term on the 



same holds for vertex- and box-like contributions. The processes appearing on the rhs of 
Eq. (1A.2I) are shown in Fig.s[22]and[23l The advantage of this formulation is the following: 
all the PT-rearranged (squared) amplitudes appearing on the rhs are at least one loop lower 
than the amplitude on the Ihs. Therefore, one can actually reconstruct the Ihs, by working 
directly on the rhs, because one knows how to pinch at lower orders. 

To see how all this analysis makes finally contact with the main question at hand, 
namely which way to pinch in higher orders, let us focus on Fig. [22l There it is clear 
that the product involves the PT-rearranged one-loop on-shell amplitude for the process 



q{pi)q{p2) 9{ki)g{k2), whose construction is abso^ 



been described in great detail in the literature 



34 



ute. 



y fixed and well defined, and has 



46| In fact, it was the first explicit 



example [46j demonstrating the universality (process-independence) of the PT gluon self- 
energy: the resulting gluon self-energy does not depend on the embedding process (quarks 
to quarks or quarks to gluons or gluons to gluons, etc). The PT-rearranged one-loop 
(l{Pi)(l{P2) —>■ g{ki)g{k2) is obtained following exactly the same PT procedure as for the 
process with only quarks as external particles. In particular, the three-gluon vertices in 
graphs (6) and (c) in Fig. [22] must be exactly as shown, i.e., the one injected with q has 
undergone the PT decomposition (and has become F^), while the ones injected with ki and 
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FIG. 23: The PT-rearranged tree-level amplitude for the process qq ggg; squaring and further 

I I 

pinching triggered by the longitudinal momenta inside the polarization tensors [30] will furnish the 
second term on the rhs of Eq. ()A.2p . 

k2 remain unchanged. Let us now return to the two representative two-loop diagrams of 
FigEH After their PT rearrangement, the two-particle Cutkosky cut on (a) denoted by (02) 
in Fig. [21] must reproduce (c) ® [(/) + (g)] in Fig. [221 and cut (04) on (6) must reproduce 
(b) (g) [(/) + ((?)]. Obviously, if we were to modify the internal three-gluon vertices of (a) 
or (b) in Fig. [21] in any way, this identification would not work: one must modify only the 
vertex injected with q (turning it to F^). This argument may be generalized to include all 
remaining two-particle and three-particle cuts, making the above conclusion completely air- 
tight. Under the light of these observations it should be clear why, for example, the relevant 
full three-gluon vertex entering into the SD equation for the gluon self-energy {viz. Figj6]) is 
indeed T^^^ (constructed in subsection II VBI ) and noi the vertex T^^^. The latter could not 
be consistently constructed inside loops, because any preferred direction (i.e. the direction 
determining the "would-be" background field) is immediately at odds with the unitarity-cut 
arguments developed above. 

The arguments presented here do not postulate at any point the existence of any relation 
between the PT and the BFM. On the other hand, on hindsight, all conclusions drawn (for 
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example, the T^^^ versus T^^^ issue) are in complete agreement with the known PT/BFM 
correspondence. Specifically, switching now to the BFM language, the fact that internal 
vertices should not be touched is precisely what the unique set of BFM Feynman rules 
dictates: since one cannot have background fields propagating inside loops, all internal 
vertices have three quantum gluons merging. This is exactly what one finds, e.g., when 



computing the two-loop gluon self-energy [l7|: as a subset of the calculation one will have 
to consider the one-loop vertex T^^\ but will never encounter the one-loop vertex F^il-- 
(constructed in [l5| and studied in |29!]). 

APPENDIX B: A BRIEF DISCUSSION OF RENORMALIZATION 

This appendix is meant to outline the general framework for dealing with the issue of 
renormalization in the context of the PT. The emphasis is put on the various conceptual and 
methodological issues involved, rather than an explicit proof of renormalizability. In par- 
ticular, we consider this discussion necessary for convincing the reader that renormalization 
poses no problem whatsoever for the self-consistent implementation of the PT. 

The analysis presented thus far assumes implicitly that the theory is renormalizable (as 
is QCD in (i = 4) or superrenormalizable (as is QCD in = 3), and that all momentum 
space integrals have been regularized by resorting to a regularization scheme that preserves 
the gauge symmetry (obviously there is little point in applying the PT to a theory that 
is ill-defined to begin with). Specifically, throughout the paper we have adopted the most 
widely used such scheme, namely dimensional regularization. Given that the original theory 
is renormalizable (by assumption) it should be clear that there is no step throughout the 
PT procedure that could jeopardize renormalizability. Indeed, all that the PT really does 
is to trigger STIs. The latter are a direct consequence of the original BRST symmetry of 
the theory; therefore, within a suitable regularization scheme (such as dimensional regular- 
ization) they will be preserved by renormalization (z.e, they will not get deformed). It is 
important to emphasize that the latter property holds true for the BQIs as well; they too 
are a consequence of the BRST symmetry and (under the same assumptions) do not get de- 
brmed either. Notice that this is completely different from the case of the Nielsen identities 



481], describing the gauge fixing parameter dependence of the bare Green's functions (we do 



not use them here). In this latter case, one needs to extend the BRST symmetry to include 
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the variation of the gauge fixing parameter. This, in turn, will spoil the original BRST 
invariance of the theory, impl ying that the latter identities get deformed by renormalization 



already at the one- loop level 43 1. 



For concreteness let us assume that we start the PT procedure in the renormalizable 
(linear) Feynman gauge (RFG), as we have done throughout this article. Let us denote by 
Za the gluon wave-function renormalization, by Z3 the vertex renormalization constant for 
the three-gluon vertex Vap,v^ by Z2 the usual ghost wave-function renormalization, and by 
Zx the ghost-gluon vertex renormalization constant. Notice also that, the BRST symmetry 
demands that Z^/Za = Z1/Z2. Then, the fundamental STIs employed when carrying out 
the PT survives renormalization, simply because all counterterms necessary to render it 
finite are already furnished by the usual counterterms of the RFG Lagrangian. This is, of 
course, a direct result of the basic assumption the the theory in the RFG is renormalizable: 
once all counterterms have been supplied in the RFG, the STI which is studied in the same 
gauge, will continue being valid. 

The actual implementation of the renormalization procedure proceeds along the lines 
described in 30|| for the two-loop case. One should start out with the counterterms that 
are necessary to renormalize individually the conventional Green's functions in the RFG. 
Then, one should show that, by simply rearranging these counterterms following the PT 
rules, one can renormalize the PT Green's functions. Notice also that, due to the validity 
of the Abelian Wis, the renormalization constants before and after the PT rearrangements 
are related to the gauge coupling renormalization as follows: 

Z2 rv2rv — 2r7—l 1y2ly — 2ly — l V — 1 I'D 1^ 

g — ^1^2 —^1^2 — ^A ■ V^-^) 

After rearranging the original RFG counterterms in such a way as to render the PT 
Green's functions finite, one should be able to verify that the resulting counterterms are 
in fact identical to those obtained when carrying out the BFM renormalization program as 
explained by Abbott [l^, i.e. by renormalizing only the background gluons, the coupling 
constant (7, and the quantum gauge-fixing parameter ^q. Thus, the relevant renormalization 
constants are given by 

g^ = Zgg, A, = zfA, eQ = Z^^^Q, Z^^ = Za. (B.2) 

The renormalization of is necessary due to the fact that the longitudinal part of the 
quantum gluon propagator is not renormalized. As pointed out by Abbott, in the context of 

65 



the BFM this step may be avoided if the calculation is carried out with an arbitrary rather 
than the BFG .^q = 1. Of course, as we have seen, the PT brings us effectively to C,q = 1; 
thus, when interpreting the resulting counterterm from the BFM point of view, one should 
keep in mind that gauge-fixing parameter renormalization is necessary. The renormalization 
of not only affects the propagator-lines, but also the longitudinal parts of the external 
vertices; it renormalizes precisely the F^ part, as can be seen from the corresponding BFM 
Feynman rule for the three-gluon vertex (see Appendix [F]) . 

All the above ingredients must be combined appropriately in order to demonstrate the 
renormalizability of the new SDK; we shall not pursue this point any further. 



APPENDIX C: FADDEEV-POPOV EQUATIONS 

As a first example of the use of the FPE introduced in Section llllDt let us differentiate the 
functional equation fl3.37p with respect to the ghost field c*; after setting the fields/anti-fields 
to zero we get (relabeling the color and Lorentz indices) 

Fc™c"(g) + iqTcmAi^^iq) = 0, (C.l) 

which can be used to relate the auxiliary function TcmAi^iq) with the full ghost propagator 
D°'''{q). Due to Lorentz invariance, we can in fact write Tcmji*n[q) = qu^c^A'^^^i^), and 
therefore 

Tcmcn{q) = -iqTc^At"{q) = -iq'^TcmA*«{q). (C.2) 

On the other hand, due to our definition of the Green's functions [see Eq. (13.181) ]. one has 
that 

zD^'-{q)T,r,4q)=6^^, (C.3) 
and therefore we get the announced relation: 

= q4q'D^-iq)]-\ (C.4) 

As a second example, let us differentiate Eq. (13.371) twice, once with respect to and 
once with respect to c^, and then set the fields/anti- fields to zero; in this way we get the 
identity 

TcM^e™ {k, q) + «g^rcM^A™ {k, q) = 0, (C.5) 
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which is particularly useful for the PT construction. Notice that this identity is the equiv- 
alent of the one introduced in Section III Al relating the conventional H function with the 
trilinear gluon-ghost vertex. All these identities can be easily checked at tree-level; for 
example, using the Feynman rules presented in Appendix [0 we have 

Differentiation of the functional (??) with respect to a BFM source f2 and a quantum 
gluon field A or a ghost field c and a background gluon A, provides instead the identities 

{ki + k + q = Q) 

TnrA^cmik,q) +iqf'TnrA^A-';r.{k,q) = gf^'^^'gup, (C.7) 
r,,.i„g™(A:, q) + 2g^r^.;f„^™(A;, q) = -ig T'^'T ,r ^^-k^) , (C.8) 

that can be easily checked at tree-level. 



APPENDIX D: SLAVNOV- TAYLOR IDENTITIES 

STIs are obtained by functional differentiation of the STI functional of Eq. fl3.30p with 
respect to suitable combinations of fields chosen following the rules discussed in Section [III CI 



1. STIs for quark proper vertices 

We begin by deriving the STI satisfied by the trilinear quark-gluon vertex (see, e 



■9: Q, 



33(1). From our general discussion of Section [III C| for obtaining this identity we need to 
consider the functional differentiation 



6c^{q)6ilj{p2)6ilj{-pi) 
from which we obtain the equation 



= q + p2=p^, (D.l) 



rc-yi7(-g)rAdV'V^(P2, -Pi) + T^,^^m{-pi,q)T^^{p2) + T^^{pi)T^^m^,{q, -pi) = 0, (D.2) 
where the two-point function r^^(p) is defined through the identity 

tS{p)T^^{p) = I. (D.3) 
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Using then the relation of Eq. (lC.4p . we get the STI in its final form, namely 



9"r^a^^(P2, -pi) = {r^,^,a'(-pi,g)r^^(p2) + r^^(pi)r^ca'^*(g, -pi)} . (d.4) 

The three-point auxiliary functions appearing in the equation above can be constructed 
using the Feynman rules reported in Appendix [Fl At tree-level the above identity can be 
trivially checked, with the rhs being proportional to two inverse quark propagators 

= ge [{^i - m) - (^2 - m)] . (D.5) 

Contrary to the case of QED, where this generalizes directly to all orders, in the QCD case 
it is valid only to lowest order, due to the non-linearity of Eq. flD.4p . 

The STI satisfied by the quadrilinear quark-gluon vertex (induced beyond tree-level) can 
be derived by considering the following functional differentiation 



= h + k2+P2=Pi. (D.6) 

$ a>*=0 



6c"'{ki)6A]^{k2)6ilj{p2)6'^{-pi 
Carrying out the functional differentiation and using again Eq. (lC.4p . we get 

ki'^A]^A^^^ik2,P2,-Pi) = [klD'^'^' {ki)]^^T,^,^^,^' {-pi, ki)T ,^n^^{p2, -P2 - k2) 
+^A^^i,{Pi - h, -pi)r^^™'^*(A;i, -P2 - ki) + r^m'^n^-7(/c2, -ki - k2)TAd^^{p2, -pi) 

+'^A^ri^cm'{p2, -Pu h)T^^{p2) + r^v^(pi)r^nvc'"'v;*'^P2, 

+rc-'A7^^(^2,P2, -Pl)rAdA;)(^2)}- (D.7) 

Clearly this identity starts at the one- loop level (recall that Fs represent IPI functions). 



2. STIs for gluon proper vertices 



Let us start by deriving the well-known STI for the trilinear gluon vertex 
considering the functional differentiation 

535(F) 
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33|. By 



6c''{q)6A]^{ki)6A]^{k2 







q + ki + k2 = 0, 



(D.8) 



$,$*=o 



and using Eq. ( ]C.4I) one obtains 

q"TAaA'pA';}iki,k2) = [g^/^""'(g)] {T^a'AnA*;<{k2,ki)TAdA^{ki) + T^a'A^A'^'<ikuk2)TAaA^{k2)^ . 

(D.9) 
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Notice that since we are working with the minimal action (see subsection IIII Cp . one has 

ri°i^.(g) = ^g'5'"'Pa/3(g), (D.io) 

and therefore 



(D.ll) 



As an overall consistency check of our definitions and conventions, notice that Eq.s (12 .Tp 
and fID.llI) imply that — r^a^b (g) = 6"'^Ila/3{q), in full agreement with Eq. (I3.18p . Using the 
above relation, we can now check the identity at tree-level; we get 



q-T^ZA^^,ik,,h) = [g^/^'^'^'(g)](°Mr(2^„^,,(A;2,A:i)r^°]^.(A:i) 

+ rfj^j,^..(A:i,fc,)ri°]^„(fc,)} 

= ^^/""" [{9f.ukl - h^h,) - ig,,kl - . (D.12) 



Notice also that Eq. ( ID. lip allows us to compare the STI of Eq. ( ]D.9p with that of Eq. (I2.10p . 
which is written in the conventional formalism; in this way we get the identity (factoring 
out the color structure) 

H^yih, /cs) = TcA^A'ih, h)- (D.13) 

We pause here to show what would have happened had we worked with the complete 
generating functional. In this case, due to the extra term appearing in the master equation 
(I3.27P satisfied by the complete action, the differentiation carried out in Eq. (ID.SP would 
generate two more terms with respect to the ones already appearing in Eq. flD.9p . namely 

5'^"A;2.r,.^™g,(A;i, h) + S'^^'h^T^aAn-Ah, h). (D.14) 

To get to the terms above we have used the equation of motion of the Nakanishi-Lautrup 
multiplier B eliminating the latter in favor of the corresponding gauge-fixing function JF. 
Then, making use of the FPE (IC.Sp . we get 

- i6'^'^k2uk2^'rc'^A"-A*."' {h, h) - iS'^"^ki^ki^TcaAnA*->{k2, ki), 



:D.15) 



so that we finally would get the STI 

2 j-)aa' I 



q^'VA^^Ar^A-Ski.h) = [q'D''^\q)]\V^.,j,^A.-,{k2,h) V'j.d A^{k{) - id''"' k^^k^^ 



dn 1 



{ki, k2) r^d^„(fc2) - i6 '^k2-^k2u 



(D.16) 
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where we have indicated exphcitly that the two-point functions are to be evaluated from 
the completed functional (for the three point functions appearing in the STI above there is 
no difference). We then see that the difference amounts to a tree- level piece appearing in 
the two-point function, as has been anticipated in our general discussion of subsection IIII CI 
(recall that we are using the Feynman gauge ^ = 1). In particular notice that we correctly 
find the relation (g) = {A-^)-'Jq). 

Another STI that will be needed in the PT construction is the one involving the quadri- 
linear gluon vertex; carrying out the functional differentiation 



5c^{h)6A^{k2)6A;{p2)SA^^{-pi 
and using Eq. (1C.4|) . we arrive at the result 



= h + k2+P2=Pu (D-17) 

$ $.=0 



^frA-A;jA-Aj(fc2,P2, -Pl) = [fci£'"""'(/Cl)]|r^^'^.^.7(-pi,fc2 +P2)rAdA^A-(fc2,P2) 
+'^cm'ArA*-'iP2,k2-Pl)'^A'iA^A%{k2,-Pl) + (/i;2, ^2 " Pi)r^d^r^. (^2, "Pi) 

+rc-'A-A»A7(P2, -Pl, k2)Ty^dAu{k2) + ^ c"-' A^A^^A'^' (^2, -Pl, P2)'r AdA-iP2) 
+rc-'AUA^A7(^2,P2, -Pl)rAdAj(Pl)}- (D-18) 

At tree-level, notice that only the first two lines of this identity are different from zero; then, 
using the Jacobi identity, we obtain 

k'ir%^,^^r,s{k2,P2,-Pi) = [fc?Z^'"'"'(A:0](°){rfi^,^.,(-pi,fc2+P2)ri°]^„^,.(A:2,P2) 

+ (^2, P2 - Pl)'^%ArAs^ (P2, -Pl) } 

= 1 {g^^k,, - g,„k^,) + r^r- {g,,k,, - g,„k,,) 

+ r'^'n (g.ph. - 9uak,,) }. (D.19) 

3. STIs for mixed quantum/background Green's functions 

Let us consider a Green's function involving background as well as quantum fields. 
Clearly, when contracting such a function with the momentum corresponding to a back- 
ground leg it will satisfy a linear WI [see, e.g., Eq.s (12.371) - (I2.40p ]. whereas when contract- 
ing it with the momentum corresponding to a quantum leg it will satisfy a non-linear STI. 
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= q + ki + k2 = 0, (D.20) 

<s>,<s>*,n=o 



Let us then study the particularly interesting case of the STI satisfied by the vertex T^^^ 
when contracted with the momentum of one of the quantum fields. Taking the functional 
differentiation 

we get 

(D.21) 

Notice that the same result can be achieved by contracting directly the BQI of Eq. flE.24p 
with the momentum of one of the quantum fields and then using the STI of Eq. flD.Op 
together with the BQIs of Eq.s (lE.2p and (]E.26P to bring the result in the above form. 

It is particularly important to correctly identify in the above identity the missing tree- 
level contributions (due to the use of the reduced functional, see also the discussion in 
Section [E2|) . In order to do that, one can either work with the complete functional and use 
the FPE fIC.Sp . or add them by hand using Eq. flE.24p . obtaining in either cases the STI 



- t9n''{klg^.-hak2u). (D.22) 



This STI can be further manipulate by using Eq. fID.lip and the FPE (IC.SP for rewriting 
the term proportional to r^A(^2) as 

r^™';ia^«(g, k2)T AiA^{,k2) = r^™'ia^«(g, A;2)(A"^)^"(A;2) + /i;2i.r^™'iag„(g, k2) 

+ «^?r^"A;2.r,w^,e(-A;i). (D.23) 

On the other hand, employing Eq. (lC.4p we find 

[klD"'"^\k^)]i^g^'^k2.)T^^,^,^^{-k^) = -ign''k,^k2:.- (D.24) 

so, inserting Eq. flD.23P back into Eq. flD.22p we see that the term above partially cancels 
the tree level contribution, thus leaving us with the STI 

k^iT^^^^^^{k,M) = [kfD"^^'ik^)] \T^^,^^^,4k2,q)T^^^^iq) + T^^,^^^^,iq,k2)iA-')t:ik2) 
+ k2,T^^, ^.^_„(g, k2)] - igr'^'^klg^,. (D.25) 
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4. STIs for the quark SD kernel 



In addition to the STIs for IPI (proper) vertices, the PT construction for SDEs requires 
the additional knowledge of the result of the action of longitudinal momenta on connected 
kernels. The first one of these kernels is encountered in the construction of the PT gluon- 
quark-quark vertex, and can be written as follows (see Fig. [TT!) 



+ i^A]i^^i^',-Pi)'>'S{i')iTArn^^{p2,-£'). (D.26) 



where i = k2 + P2 = Pi — ki and = /ci + P2 = Pi — ^2- Then, using the STI of Eq. flD.4p 
and the relation (ID. 30 . one gets the results 

k^,iT^^^^{i,-p,)tS{i)tT^r.^^ip2,-i) = -[klD^^^'ik^)] {r^,^,w(-pi,fci) 

+ ^r^^(pi)r^c'"'v;*(^i' -Pi)'S'(^)} r^j^^(p2, -^), 

(D.27) 

k^tT^r.^^{e',-p,)zS{e')tT^r.^^{p2,-i') = -T^.^^{e',-p,)[klD"^"^'{k,)] {T^^^,^,{k,,-n 

+ tSii')T^.^^^,i-i',k,)T^4p2)}. (D.28) 

We then see that the first term in Eq.s (ID. 2711 and (ID. 2811 will cancel the first two terms of 
the STI of the IPI vertex of Eq. (1D.7|1 . and we finally arrive at the STI 

ki^A]^A^^^{k2,P2, -Pi) = [fci (fci)] |r,„.'^n^*7(A;2, -ki - k2)Tj^d^^{p2, -pi) 

+ r^v;(Pi)^A^vc'"'Vi*(^2, ki, -pi) + ^A^^*V^c™'(P2, -Pi, ki)T^^{p2) 
+ ^cm'A';'^i,{k2,P2,-pi)TAdA^{k2)^ , (D.29) 

where we have defined the auxiliary kernels 

^A^^cm'^*{P2,ku-Pi) = r^„^,™'^*(p2,A;i,-pi) 

+ iT^^m'^*{ki, -pi)iS{i)ir^„^^{p2, -i), (D.30) 

^A?v*v;c'"'(P2,-pi,fci) = r^^v*v;c'"'(P2,-pi,^i) 

+ zr^„^^(£,-pi)z5'(£)zr^.^-,„,(-£,fci). (D.31) 
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5. STIs for the gluon SD kernel 

In the construction of the SDEs for the gluon self-energy and three-gluon vertex, one 
needs the knowledge of the STI satisfied by the kernel (see Fig. 

+ i^A%AY}A<i {h , (£)zr^e'^„^. (^2, P2) 

+ iVAs^A«^At{k2,t)i/\i,{(i')iV^.,^^^^{h,p2). (D.32) 

Using the above relation, together with STI of Eq. (lD.9p . we find the following result 

KlTAs^A^At{h.e)l/\i,{i)lT^,j^^^j^^{k2,P2) = -[klD^'^'ik,)] X 



-pi,£)P"'(£)+2r 



{£, -pi)TAdAs„iPl)^Te'ii) \T^,, ,„^r{k2,P2) 



(D.33) 



In this case this is, however, not the end of the story, since the first term in the equation 
above still contains (virtual) longitudinal momenta, which will trigger the STI of Eq. flD.Op 
together with the FPE (IC.Sp . After taking this into account, we obtain 



A:r^r^jA^Ae(fci,£)2A«,(£)2r^e;^„^.(A;2,p2) = -[A;fD™™ {k,)] x 

x{ ^c-^'A%Alfi-Pl^^)+^'^crn'AtA';<{^^-Pl)^Ad^A%{Pl)Ki{^) ^ A-J A^ Arik2 , P2) 
+iT^^'As^A-Pu^)D^^\^) r^^'A-A^fe, k2)TAd^Ai;{h) + r^.'^„^*7(A;2,P2)rAdA^(P2) 



(D.34) 



Similarly we find 

k'^zTAs^A^At{h,i')iAli,{i')zT^.,^^^^{k,,p2) = -[A;2d-™'(A;i)] x 



y<\TAiA^At{k2,n T^^,^.,,Jp„-i')+tAZ,{i')T 



c^'a^Ia*/^ 



?',P2)rAdA-(P2) 



+iD^^\i') T^,A^A*Ak2, -pi)TAdAs{pi) + T^eAsA*-r{-pi, /Ca)!^^^!? (^2) 



(D.35) 



As before, after combining these results with the four-gluon IPI vertex STI of Eq. (ID.lSp 
we arrive at the needed STI for the four-gluon SD kernel, namely 

ki}CA^A])ArA%{k2,P2, -Pi) = [fci (fci)] I r^„/AnA.7(/c2, -ki - fca)!^^^-^^ (P2, -Pi) 
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+ ^c^'ArA%A';<{P2, -Pl, h)'^A^A^{k2)y 

where the following auxiliary kernels have been defined 

u a d 1^ <^ a 

+ i^A%A^Ai {h, (^')^r^„/^e'^*7 (-^', P2) 

,7(^2,^2,-^1) = r^w^„^.^.7(A;2,P2, -Pi) 



/c 



'(^,-Pl)«^ee''(^)^r^<A'-A"(^2,P2) 



/c 



7(^2,-^1,^2) 



rc-'A^AjA7(P2,-Pl,/i;2) 



+ «rgm'^sge'(-pi,^)^i^'''' (£)irg.'^.^.7(p2, ^2) 

c Pa 



(D.36) 



(D.37) 



(P2,-^'), (D-38) 



(D.39) 



APPENDIX E: BACKGROUND-QUANTUM IDENTITIES 



BQIs are obtained by functional differentiation of the STI functional of Eq. fl3.32p with 
respect to combinations of background fields, quantum fields and background sources. 



1. BQIs for two-point functions 



The first BQI we can construct is the one relating the conventional with the BFM gluon 
self-energies. To this end, consider the following functional differentiation 

6^s' (r') 



Sn^^{p)5A^f,iq) 

6^s' (r') 







which will give the relations 



•s>,<s>*,n=o 



iT 



Ab 



(?) 



^gls^' + rn^^A^-^iq) 



q + p = 0, 
q + p = 0, 



^A^Aijiq)^ 

^Ad^Al^iq)- 



(E.l) 
(E.2) 

(E.3) 
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FIG. 24: Expansions of the gluon anti-field and BFM source in terms of the corresponding com- 
posite operators. Notice that if the anti-field or the BFM sources are attached to a IPI vertex, as 
shown in the first line, such an expansion will in general convert the IPX vertex into a (connected) 
SD kernel. The equivalence shown is therefore not valid at tree-level {e.g., in the case of three-point 
functions such an equivalence would imply that the kernels shown on the rhs of the corresponding 
expansions would be disconnected); when present, the tree-level needs to be added by hand, as 
explicitly shown in the two expansions of the second line and the last one of the third line. This 
type of expansion allows one to express the terms appearing in the BQIs in a form that reveals 
kernels appearing in the STIs [see, e.g., Eq.s (|ET4]) . (jKTS]) . (jElS]) and (jET9l) ] 



We can now combine Eq.s (1E.2P and (1E.3|) such that the two-point function mixing 



background and quantum fields drops out, to get the BQI 

- ^rf^a^.7(g)r^d^e(g)rf^6^„(g) 

= ^r^gA^(9) + 2rf^.^.7(g)r^d^6 (g) - «rf^a^-^(g)r^d^.(g)rf^.^*.(g), (E.4) 

where the last identity is due to the transversality of the Taa two-point function. 

In order for our PT procedure to be self-contained, it is important to express the IPI 
auxiliary Green's function involved in the various STIs and the BQIs in terms of kernels 
that also appear in the relevant STIs. The key observation that makes this possible is that 
one may always replace an anti-field or BFM source with its corresponding BRST composite 
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An , AS An AS 




FIG. 25: Skeleton expansion of the kernel appearing in the SDE for the auxiliary function TcAA* 



operator. Thus, for example, one has (see Fig. [2 



„r(o) 

c-'A"',AV 



,r(o) 



(E.5) 
(E.6) 



where ki and ^2 are related through k2 = q — ki. In this way we get the following SDEs (see 
again Fig. 



-F. 



(q) = -s'-'q.-Kr^Aviq) 



ik,q) 



'i^n%Aic'^{k,p) 



-6^^q^ + gf-'-'g2, [ D^'^{k,)^yl{k,)T ^^A-^AhM), (E.7) 
igr''-gZ-igr'''''gl, I D^^'(fci)A:;'(A;2)/C,.^n^j,-e(fc, fc^, ^i), (E.8) 
gr'^'gcw I D'^\k,)A^;Uk2)T,.A^.A',Ah,-q), (E.9) 

J ki 

-igr'^9.,-igr'^'go.u' [ D^\k^)A^;''^{k,)]C,.^r.^d-Ak2,k,p). 

J ki 



(E.IO) 



The kernel JCcaas appearing in the SDEs flE.81) and flE.lOp is shown in Fig. [25] and reads 



^c''A^Ajc'=(^5 ^2; ^l) — ^c'^A]iA^c'^{k, k2, ki) 

+ i^Ai;A-^A-Ak2, -k - k2)iA'^'^,{k + A;2)iF^m^r'ge(A; + k2, ki) 

p' 

+ iTcaAi^c^{k2,-ki- k2)iD'"^'{ki + k2)iT^r'j^ni:e{k,ki). (E.ll) 

Finally notice that the auxiliary function F^^^* corresponds precisely to the auxiliary func- 
tion Aq,/3 introduced in Eq. (12.331) . and therefore its part proportional to gaf3 corresponds to 
the scalar function G{(f). 
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2. BQIs for three-point functions 



The relation between the trihnear quantum gluon-quark vertex and the trihnear back- 
ground gluon-quark vertex, can be obtained by considering the following functional differ- 
entiation 

^ ' -0 g + p2=Pi. (E.12) 



We then get 



+ r^*^ng(-pi,9)r^^(p2) + r^v;(pi)r^ng^*(?>-pi)- (E-13) 

In order to explore further the all-order structure of these two auxiliary Green's functions, 
replace the BFM source with the corresponding composite operator using Eq. (]E.6[) . thus 
obtaining 

iT^ni^,{q,-pi) = ^r[^°]^„;g,„/ ^ zL'™'"'(fci)zA^;;'^(fc2)/C^n^c-v;* (P2, ^i, -pi), (E.14) 

i'^^'M-i-Puq) = ^r[°^.n'-„' / ^^"''™(^l)^A^^''n(^2)/C^n^,^c-(P2,-Pl,^l)- (E.15) 



= g+p2=Pi, (E.16) 

<s>,<s>*,n=o 



where the kernels /C^n^^™^* and }Cj^„^*^^m have been defined in Eq. flD.30p and flD.3ip . As it 
is clear from the two equations above, while the (auxiliary) functions appearing in the STIs 
and BQIs ought to be IPI (Ihs of the equations), that is not true for the kernels appearing 
after using the substitutions of Eq. flE.5l) or flE.6l) . which, in fact, consist of both IPI and 
IPR diagrams (rhs of the equations). 

For the BQI involving the three-gluon vertex a similar result can be obtained; choosing 

6^S'{T') 
SQ%{q)6A;{p2)6A^A-P 
we will get 

+ ^n-A^^A*;'{-PuP2)^Ai^Ar{P2) + rngA-A7(P2, -Pl)'^ A^A^^iPl) ■ 

(E.17) 

Again we can write 

^rnaAjA7(-Pl,P2) = / tD^'^ih)tA^2ik2)JC,-A^AsA:Ak2,-pi,P2), (E.18) 

™ J ki 

^^niArAT{P2,-pi) = / ^/^"^''"(/^i)^A:;::;(/.2)/C,™^.^.^*.(fc2,P2,-pi), (E.19) 
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where the kernels appearing in the above equations have been defined in Eq.s f lD.37p 



and flDlSSj) . Notice the emergence of the pattern exploited in the application of the PT 



to the SDEs of QCD: namely that the auxiliary functions appearing in the BQI satisfied 
by a particular Green's function can be written in terms of kernels appearing in the STIs 
triggered when the PT procedure is applied to that same Green's function. The BQI of 
Eq. (]E.17p gives at tree-level the result 

^A-A-As (P2, -Pl) = rAgA-Aj(P2, -Pl)- (E.20) 

This is once again due to the use of the reduced functional: in fact in such case the two 
(tree-level) vertices need to coincide, since the difference between them is proportional to 
the inverse of the gauge fixing parameter (see Appendix [F]) and therefore entirely due to the 
gauge fixing Lagrangian. To restore the correct tree-level terms one would have to use the 
complete functional; in that case the differentiation of Eq. ( 1E.16|) shows the two additional 
terms 

- 5'^''pi„V^a^AiAP2, -Pi) + 5'^''P2p^n-A%A-PuP2), (E.21) 
which, with the help of Eq. flC.7|) become 



-i5'^'PlaPl^T^aArA*^i {P2, -Pl)-i5'^''P2pP2y^n-As^A*;' (-Pl, P2) +9r'''iqapPla+gaaP2p)- (E.22) 

Therefore we get the final identity 

^r|„^,^, (P2, -Pl) = IWIS'"'^ + rQa^*7(-g)]r^dA.^j(p2, -Pi) + gr''{qapPla + gaaP2p) 
+ ^n%A%A*^^{-Pl.P2) r^dA-fe) -«'5'^>2pP27 

+ ^n-ArA*AP2,-pi) T%j^s{pi)-i6'^'piaPi^ , (E.23) 



a p a 



which gives the expected tree- level result. Once again we see that the difference between 
working with the reduced and complete functional lies in some constant (tree-level) terms 
that one recovers after applying the FPE for writing the STI/BQI at hand in the same form 
using F or Fq. Thus, opting for the fast way of deriving the STI/BQI with the reduced 
functional and adding the correct tree-level term, we write the BQI in its final form 

^riaA-Aj(P2, -Pl) = [iglS^'^ + T^aA*',i-(l)]^Ad^A-Ai{P2,-Pl) 

+ ^Q-As^A*;'{-PuP2)TAdAr{P2) + ^ Q^ArA*;' {P2, -pOF^d^. (pi) 

+ gr'{P2pga.+Piagap). (E.24) 
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3. BQI for the ghost-gluon trilinear vertex 



In this section we are going to derive the BQIs relating the ghost sector with the BFM 
ones. We start from the trihnear ghost-gluon coupling, for which we choose the following 
functional differentiation 



= ki + k2 = q, (E.25) 

<I>,4>*,Q=0 



5(]^(-g)5c™(A;i)(5c"(A;2 
thus getting the result 

~ ^c"^A*'^{'~kl)TQay^dcn{kl,k2) —TQa^nic*d{kl,k2)Tcdcn{k2). (E.26) 

At tree-level we then correctly recover the result 

= -gn^'ih - k2)a, (E.27) 

(in this case there is no difference between using the complete or reduced functional). 



APPENDIX F: FEYNMAN RULES 
1. and BFM gauges 

The Feynman rules for QCD in gauges are given in Fig. [261 In the case of the BFM 
gauge, since the gauge fixing Lagrangian is quadratic in the quantum fields, apart from 
vertices involving ghost fields only vertices containing exactly two quantum fields might 
differ from the conventional ones. Thus, the vertices F^^^ and F^^^^ have to lowest order 
the same expression as the corresponding R^ ones F^^^ and Taaaa (to higher order their 
relation is described by the corresponding BQIs). 



2. Anti-fields 

The couplings of anti-fields with fields is entirely encoded in the BRST Lagrangian of 
Eq. (13.201) . When choosing the BFM gauge the additional coupling g f"-"^^ A*^' A^c"' will arise 
in the BRST Lagrangian £brst as a consequence of the BFM splitting A A + A. One 
then gets the Feynman rules given in Fig. [281 
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FIG. 26: Feynman rules for QCD in the gauges. The first two columns show the lowest order 
Feynman diagrams and rule respectively, while the last one shows the corresponding all-order 
Green's function according to the conventions of Eq. (j3.18p . 

3. BFM sources 

The couplings of BFM sources with fields can be derived from the Faddeev-Popov 
ghost Lagrangian, since making use of the extended BRST transformation of Eq. (13.311) we 

get 

/:PPG = -c^sT^^u 3 -c'^^7/'^"^"(3A;r)A^: = -gn-c^n^A^^. (F.l) 
The corresponding Feynman rule is finally given in Fig. [211 
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FIG. 27: Feynman rules for QCD in the BFM gauge. We only include those rules which are different 
from the ones to lowest order. As usual, gray circle on a gluon line indicates a background field. 
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